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ABSTRACT 


A  recent,  anticipatedly  exact,  analjrtical  theory  of  homogeneous 

turbulence  is  specialized  to  stationary  isotropic  turbulence  of  very  high 

Reynolds  number,  maintained  by  idealized,  low  wavenumber  isotropic  driving 

forces.  The  characteristic  wavenumber  k  =  e/v    and  Reynolds  number 

o  o 

R    =  V  k  "  /v  ,  where    r       is  the  rms  velocity  in  any  direction,     e     is  the 
ooo  o  ^  -V  J 

power  dissipated  per  unit  mass,   and     v     is  the  kinematic  viscosity,  are 

1/3 

introduced.  For  R    »  1  it  is  found  that  the  inertial  and  dissipation 

o  '^ 

ranges  extend  over  wavenumbers  satisfying  k  «  k«  R  k  .  In  this  case,  the 

ooo 

time  correlation  function  and  ensemble-averaged  impulse  response  function  for 
an  inertial  or  dissipation  range  Fourier  mode  are  found  to  be  J-,(2v  kr)/(v  kr) 
and  exp(-vk  r)j-(2v  kT)/(v  kx),  respectively.     The  energy  spectrum  in  these 
ranges  is  determined  by  a  nonlinear  integral  equation,   involving  the  time 
correlation  and  response  functions,  which  is  suitable  for  solution  by  iteration. 
The  solution  is  of  the  form  E(k)/(r^v  =  (kA^j)"-^     f(kA^),  where  k^  =  R^^^^k^ 
and  f  is  a  universal  function.     In  the  inertial  range  it  reduces  asymptotically 
to  E(k)  =  f(0)(cv   )   '  K     '    .     The  parameter  f(0)   can  be  obtained  by  quadratures, 
without  solving  the  integral  equation  for  E(k),     Spectral  energy  transport 
throughout  the  inertial  and  dissipation  ranges,  and  also  up  through  wavenumbers 
k  2:  RqJ^q*  is  found  to  proceed  by  a  cascade  process  essentially  local  in  wave- 
number  space  J   the  direct  power  delivered  by  all  modes  below  k  to  all  modes  above 
k'»  k  is  of  order  e(kA')         if  k  and  k'  both  lie  within  the  inertial  range. 

The  spectrum  contains  an  exponential  factor  for  k»k  ,  and  mean-square  velocity 

1/3 

derivatives  of  all  orders  are  finite.     For  R     '-^  »1  the  skewness  factor  of  the 

0 

distribution  of  the  nth-order  longitudinal  velocity  derivative  approaches  the 

asymptotic  form  A  R  "  '^    ,  where  A     is  a  universal  constant, 
no'  n 
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The  k     '      inertial  range  law,  and  not  the  Kolmogorov  k  law,  is 

found  to  be  followed  closely  by  the  one-dimensional  spectirum  of  the  total 

kinetic  energy  measured  by  Laufer  on  the  axis  of  fully  turbulent  pipe  flow  at 

Reynolds  number  R     '•'--'  15,     The  nature  of  the  isotropic  solution  described 

above  suggests  that  the  solution  for  anisotropic  turbulence  may  yield  anisotropic 

asymptotic  partition  of  energy  for  wavenumbers  in  the  inertial  and  dissipation 

ranges  because  of  anisotropy  of  the  time  correlation  and  response  functions  for 

k  <;    R  k  .     This  offers  a  possibility  of  explaining  the  high  wavenximber  anisotropy 

found  by  Laufer.  It  is  found  that  the  asymptotic  dissipation  range  scaling  given 

by  the  present  theory  is  better  supported  by  grid  turbulence  results  of  Stewart 

1/3 

and  Townsend  (R  '•^~  2,3-3.7)  than  the  Kolmogorov  scaling, 
o 

It  is  stressed  that  extreme  caution  must  be  exercised  in  interpreting 
the  apparently  favorable  experimental  evidence  as  support  for  the  theory. 

A  comparison  is  given  between  the  present  theory,  the  Kolmogorov  theory, 
Heisenberg's  eddy-viscosity  theory,  and  the  analytical  theories  of  Heisenberg 
and  Chandrasekhar.  The  origin  of  the  disagreement  with  Kolmogorov 's  conclusions 
is  that  on  the  present  theory  the  energy  containing  modes  do  not  have  the  trivial 
convective  action  of  a  uniform  velocity  field  on  the  high  k  dynamics.  Instead, 
they  mix  weakly  statistically  interdependent  high  k  modes  into  each  other, 
thereby  importantly  affecting  the  triple  phase  correlations  which  determine  the 
mean  energy  transport  within  the  high  k  region.  The  theories  of  Heisenberg  and 
Chandrasekhar  based  on  assuming  zero  fourth-order  cumulants  may  be  obtained  from 
the  present  theory  by  discarding  some  of  the  dynamical  coupling  terms  in  the 
statistical  equations  of  motion  and  ignoring  the  effect  of  the  dynamical  coupling 
on  the  laodal  response  functions.  These  theories  also  are  inconsistent  with  the 
Kolmogorov  theory.  A  modification  of  Heisenberg's  heuristic  theory  is  examined 


-  lii  - 

which  involves  a  rational  expression  for  the  eddy  viscosity  and  leads  to  a 

-3/2 

k  '  inertial  range  law.  The  energy  transport  expression  given  by  the  eddy- 
viscosity  theory  (modified  or  not)  is  found  in  several  respects  to  be  less 
plausible  physically  than  that  of  the  present  theory.  A  brief  comparison  is 
given  with  the  theories  of  Proudman  and  Reid,  and  Tatsumi,  It  is  concluded 
that  these  theories  do  not  yield  an  inertial  range. 
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1.  INTRODUCTION 

A  theory  recently  has  been  proposed  (Kraichnan,  1958,  hereinafter 
cited  as  Paper  I)  which  is  anticipated  to  provide  an  exact  and  (in  principle, 
at  least,)  complete  statistical  description  of  turbulence  that  is  homogeneous 
in  at  least  one  direction.  The  theory  was  explicitly  worked  out  in  Paper  I  only 
for  the  analytically  simplest  case  of  stationary,  completely  homogeneous  tur- 
bulence, for  which  stationary  homogeneous  driving  forces  must  be  invoked.  In 
the  present  paper,  a  further  specialization  is  made  to  stationary  isotropic 
turbulence.  The  theory  is  employed  to  yield  the  spectral  structure,  in  both 
frequency  and  wavenuiriber,  of  the  inertial  and  dissipation  ranges  for  very  high 
Reynolds  numbers.  These  asymptotic  results  are  compared  with  experiment,  and 
the  theory  itself  is  examined  in  relation  to  current  theories.  An  attempt  is 
made  to  bring  out  as  fully  as  is  practicable  the  physical  significance  of  the 
basic  hypotheses,  analytical  structure,  and  predictions  of  the  theory.  Reference 
should  be  made  to  Paper  I  for  the  fundamental  statistical  mechanical  aspects  of 
the  theory  and  for  the  derivation  of  the  basic  statistical  equations  of  motion 
which  form  the  starting  point  of  the  present  investigation, 

2.  FUNDAMENTALS  OF  THE  THEORI 

2,1,  The  weak  dependence  hypothesis 

We  may  start  by  considering  a  uniform  incompressible  fluid  confined  in 
a  cubical  box  of  side  L.  The  velocity  field  u.(x,t)  within  the  box  may  be 
expanded  in  the  usual  fashion  in  a  spatial  Fourier  series 

ik-: 


(2.1)  u.(x,t)  =  2:e^'V(k,t), 

1  "       k        i 
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and  the  coefficients  u(k,t)  may  be  considei-ed  the  dyn?.mical  degrees  of  freedom 
of  the  fluid.     In  terms  of  them,   the  Navier-Stokes  equation  may  be  written 

(2.2)  u.(k,t)  +  vk%(k,t)  =  -  ik  P,,(k)     H      u,(p,tk  (q,t), 
where  the  dot  indicates  time  differentiation  and  the  operator 

P.  .(k)   =  6.  .  -  k'\.k. 

expresses   the  action  of  the  pressure  forces  in  maintaining  the  incompressibility 
property 

(2.3)  k^u^(k,t)  =  0. 

The  kinematic  viscosity  is  denoted  by  v.  The  summation  in  (2,2)  nins  over  all 
wavevectors  p,  q  allowed  by  the  boundary  conditions,  so  that  each  Fourier  mode 
is  coupled  to  essentially  every  other  mode. 

Let  us  assume  a  state  of  turbulent  motion  such  that  away  from  the  walls 
the  flow  is  statistically  homogeneous  in  the  sense  that  properties  defined  as 
averages  over  sufficiently  large  volumes  do  not  vary  with  location  of  the  volumes 
in  the  box.  This  implies  that  L  must  be  much  larger  than  characteristic  length 
scales  of  the  turbulence j  for  exact  homogeneity,  the  effects  of  the  walls  must 
be  vanishingly  small,  and  we  must  have  L  infinite.  Let  us  then  imagine  a  sequence 
of  boxes  of  increasing  L.  As  we  go  to  larger  and  larger  members  of  the  sequence, 
the  spacing  between  adjacent  wavevectors  allowed  by  the  boundary  conditions 
decreases;  the  density  of  allowed  modes  in  k  space  is  "^  L*^,  For  L  infinitely 
large  we  may  expect  that  an  infinitely  large  number  of  u(k,t)  are  excited  in  any 
given  wavevector  range,  no  matter  how  small.  This  is  true  regardless  of  how 
small  may  be  the  Reynolds  number  of  the  turbulence  based  on  the  rms  turbulent 
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velocity  and  the  turbulent  macroscale,  (We  may  keep  these  quantities  constant 
as  L->'^), 

As  L-><^  ,  each  u(k,t)  is  determined  by  contributions  from  more  and 
more  individual  x  space  volxiraes  of  any  given  size.  Interference  effects  of  all 
the  individual  volumes  may  be  expected  to  give  rise  to  increasingly  rapid  and 
random  oscillations  of  phase  and  amplitude  as  a  function  of  k.  On  the  other 
hand,  since  any  neighborhood  in  k  space,  no  matter  how  small,  contains  an 
infinite  number  of  excited  degrees  of  freedom,  as  L-><=^  ,  we  may  expect  that 
averages  over  all  the  u(k,t)  in  any  such  small  neighborhoods  will  yield  smooth 
results  independent  of  the  precise  sizes  or  shapes  of  the  neighborhoods.  In 
fact,  averaging  over  tiny  neighborhoods  in  k  space  is  equivalent  to  averaging 
over  large  regions  in  x  space.  The  averaging  over  tiny  neighborhoods  may  be 
replaced  by  averaging  over  a  suitable  ensemble  of  flows.  This  latter  kind  of 
statistical  average  will  be  understood  hereinafter}  it  leads  more  elegantly  to 
identical  results. 

As  we  take  the  limit  L-^  "^^  ,  the  number  of  bilinear  terms  on  the  right 
side  of  (2.2)  increases  without  limit,  but  any  given  wavevector  appears  in  at 
most  two  terms  in  the  sum  (once  as  p  and  once  as  q).  Therefore  it  would  appear 
that,  in  the  limit,  the  dynamical  coupling  among  any  finite  number  of  modes  is 
infinitely  weak,  since,  in  the  equation  of  motion  of  each  member  of  the  group, 
the  others  appear  in  only  a  finite  number  of  the  infinitely  many  bilinear  terms 
which  contribute. 

It  seems  natursil  to  assume  that  turbulence  should  be  described  by  a 
statistical  distribution  of  the  velocity  field  which  is  as  random,  in  some 
suitable  sense,  as  the  equations  of  motion  and  the  constraints  on  the  system 
permit.  In  terms  of  the  k  space  representation,  we  may  express  this  assumption 
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as  requiring  that  the  various  u,(k,t)  be  as  statistically  independent  as  is 
permitted  by  the  equations  of  motion  and  constraints.  The  above  observation 
about  the  weakness  of  dynamical  coupling  within  finite  groups  of  modes  then 
suggests  the  following  e:q)licit  hypothesis  about  the  statistical  stmcture  of 
homogeneous  turbulence: 

All  statistical  moments  formed  from  any  finite  number  of  the  Fourier 
modes  have  normalized  values  which  tend  in  the  limit  L->  °<=>  to^  values 
corresponding  to  some  distribution  in  which  the  Fourier  modes  are  completely 
statistically  independent,  subject  to  (2,3)  and  the  reality  condition  u  (k,t)  = 
■  u(-k,t).  We  shall  call  this  the  weak  dependence  hypothesis.  It  is  important 
for  the  consistency  of  the  hypothesis  that  all  other  modes  are  paired  differently 
in  the  bilinear  terms  in  the  equations  of  motion  (2,2)  for  different  given  modes. 
If  this  were  not  so,  strong  statistical  dependencies  might  arise  despite  the 
weak  d3mamical  coupling  of  the  given  modes. 

It  should  be  noted  carefully  that  weak  dependence,  which  can  be  defined 
only  in  terms  of  a  limiting  process,  is  very  different  from  independence  of  the 
Fourier  modes,  certain  properties  of  which  have  been  assumed  by  Millionshtchikov 
(I9UI),  Heisenberg  (19U8),  Proudman  and  Reid  (19$Ii),  Chandrasekhar  (1955), 
Tatsumi  (195?),  and  other  authors.  It  is  known  experimentally  that  the  joint 
probability  distribution  of  the  velocity  field  at  more  than  one  point  is  not 
normal  (Batchelor,  19^3).  In  particular,  the  skewness  and  flatness  factors  of 
the  two-point  distribution  do  not  have  normal  values.  This  implies  that  certain 
triple  niDments  of  the  u(k,t)  cannot  vanish  as  they  would  in  an  independent 
distribution  and  that  certain  fourth-order  moments  are  not  related  to  second- 
order  moments  as  in  an  independent  distribution.  As  a  complement  to  the  experi- 
mental evidence,  it  readily  is  seen  fiDm  (2,2)  that  non-zero  mean  energy  transfer 
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among  the  u(k,t)  requires  certain  non-vanishing  triple  moments,  and  it  has  been 
shown  also  that  the  equations  of  motion  imply  non-independent  relationship  of 
fourth-order  and  second-order  moments  (Kraichnan,  1957).  However,  there  is  no 
conflict  between  the  experimental  evidence  and  the  hypothesis  of  weak  dependence. 
We  may  illustrate  this  by  examining  the  skewness  factor 

<(du^/ax^)^>/<(9u^/dx^)  )'^'    ,     It   <(Su^/3x^)^>  is  expanded  as  a  sujnmation  over 
triple  moments  of  the  Fourier  modes,  it  readily  is  seen  that  the  nuniber  of  terms 
in  the  sum  which  involve  only  modes  with  wavenunibers  below  ar^y  arbitrary  value 
increases  as  L  as  L-^  °^  .  The  fact  that  the  individual  triple  moments  approach 
zero  as  the  limit  is  taken  by  no  means  implies  that  the  sum  approaches  zero.  In 
a  similar  fashion  it  may  be  seen  that  the  property  of  weak  dependence  does  not 
imply  normal  flatness  factor  for  the  distribution  of  6\L/3x,,  or  normal  relations 
among  other  moments  of  the  x  space  distribution.  Because  of  the  mutual  inter- 
ference effects  from  all  parts  of  the  box,  mentioned  above,  weak  dependence  is 
consistent  even  with  the  sharply  defined  x  space  structures  (shear  fronts  and 
vortex  filaments)  which  make  up  very  high  Reynolds  number  homogeneous  turbulence. 

2,2.  The  direct  interaction  hypothesis 

It  is  clear  from  (2.2)  that  the  total  nonlinear  interaction  of  the 
Fourier  modes  is  the  resultant  of  elementary  interactions  each  involving  a  triad 
of  modes.  As  L->-='0,  the  number  of  contributing  triad  interactions  on  the  right 
of  (2.2)  increases  as  L  .  The  hypothesis  of  weak  dependence  suggests  an  approach 
to  the  problem  of  turbulence  dynamics  in  which  each  triad  coupling  is  treated 
effectively  as  an  infinitesimal  perturbation  on  the  total  motion  due  to  all  the 
triad  couplings,  and  the  viscous  forces,  acting  together.  Such  an  approach  can 
be  carried  out  in  an  apparently  self-consistent  manner  and  used  to  evaluate  the 
physically  important  triple  correlations  of  the  u(k,t),  provided  a  second,  related 
statistical  hypothesis  is  made. 
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Let  k,  p,  q  be  three  wavevectors  such  that  k  -  p  -  q  -  0,  The  dynamical 
interaction  of  u(±k),  u(±p),  u(±q)  can  take  place  along  a  direct  path  and  along 
indirect  paths.  By  the  direct  interaction  we  mean  that  given  by  the  terms 
bilinear  in  components  of  u(±p),  u(±q)  in  the  equation  of  notion  (2,2)  for  u(±k), 
and  the  corresponding  terms  in  the  equations  of  motion  of  u(ip)  and  u(±q)j  that 
is,  the  bilinear  terms  in  the  equations  of  motion  of  these  amplitudes  which  do 
not  involve  amplitudes  associated  with  any  wavevectors  other  than  ±k,  ±p,  +q. 
By  the  indirect  interaction  we  mean  the  dynamical  coupling  among  these  three 
modes  which  would  remain  if  the  terras  representing  their  direct  interaction  were 
removed  from  the  equations  of  motion.  The  indirect  coupling  exists  by  virtue  of 
the  coupling  of  each  of  the  given  three  modes,  through  the  remaining  bilinear 
terms,  to  all  the  other  modes  in  the  system. 

Let  us  now  consider  the  triple  correlation  (u.  (k,t)u.(p,t')u  (q,t'0) • 
It  is  clear  that  the  value  of  this  correlation  is  affected  by  the  presence  of  the 
direct  interaction,  for  this  interaction  produces  increments  in  u.  (k,t)  which 
bear  a  definite  phase  relation  to  products  of  the  form  u.(p,t)u  (q,t).  We  shall 

J  — '      ax 

make  the  hypothesis  that  as  the  limit  L->  °°   is  approached  the  triple  correlation 
becomes  determined  by  the  direct  interaction  in  the  sense  that  if  the  direct 
interaction  terms  were  not  present  the  value  of  the  triple  correlation  would  be 
small  without  limit  compared  to  its  value  with  these  terms  present.  We  shall 
call  this  the  direct  interaction  hypothesis. 

Several  arguments  may  be  presented  in  favor  of  this  hypothesis,  and 
there  exists  analytical  support  of  an  a  posteriori  nature,  as  discussed  in 
Paper  I.  A  rigorous  verification  has  not  been  obtained.  For  present  purposes 
we  shall  simply  assert  the  plausibility  of  the  assumption  that  the  coupling  of 
three  degrees  of  freedom  by  the  indirect  interaction  alone  fails  to  convey  any 
phase  information  among  them  in  the  limit  where  an  infinite  number  of  weakly 
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statistically  dependent  modes  constitute  the  system-as-a-whole  through  which  the 
indirect  interaction  takes  place.     This  seems  intuitively  reasonable  if  one 
views   turbulence  as  a  mixing  or  randomizing  process.     In  any  event,  we  do  not 
present  the  hypothesis  as  an  approximation.     We  consider  both  the  weak  dependence 
and  the  direct  interaction  hypotheses  to  express  exact  properties  of  exactly 
homogeneous  turbulence  (L->°°).     There  does  not  appear  to  be  any  contrary 
evidence.     The  extension  of  the  direct  interaction  hypothesis  to  certain  fourth- 
order  moments  and  to  third-order  moments  involving  possible  external  forces 
driving  the  turbulence  is  given  in  Paper  I. 

2,3.     Perturbation  theory  determination  of  triple  moments 

We  may  evaluate  triple  moments  of  the  form  (u.    (k,t)u.(p,t' )u  (q>t")^ 
on  the  basis  of  the  direct  interaction  hypothesis  by  first  considering  a  hypo- 
thetical system  which  differs  from  the  actual  system  in  that  the  several  terms 
giving  the  direct  interaction  are  missing.     For  this  system  we  assume  the  triple 
moment  is  negligible.     Then  we  obtain  the  value  of  the  triple  moment  in  the 
actual  system  by  switching  on  the  direct  interaction  terras  as  a  perturbation. 
It  is  clear  from  what  has  been  said  that  these  few  terms  represent  an  infini- 
tesimal perturbation  for  L-><^    and  that,   consequently,   first-order  perturbation 
theory  may  be  used  without  error.     Thus,  in  the  limit  1^=^    we  have 

(2.U)  <u^*(k,t)u.(p,t')Ujj^(q,t"))   -  <5u^*(k,t)uj(p,t')tijq,t")> 

+<Ui*(k,t)6u^(p,t')uJq,t")>   +  <u^*(k,t)u^(p,t')6ujqi,t")>  , 

where  6u.(k,t)  is  the  infinitesimal  perturbation  induced  in  u. (k,t). 

In  general,  the  linear  response  to  an  arbitrary  solenoidal  infinitesimal 
perturbation  term  ?.(k,t)  on  the  right  side  of  (2.2),  switched  on  at  t  -  0,  is 
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(2.?)  6u  (k,t)  »j    ^..(kjt,3)C.(k,3)ds, 

where  ^. .(kjt,s)  is  the  infinitesimal  impulse  response  tensor.     Since  the  system 
is  nonlinear,  S^.  .(kjt^s)  is  some  (so  far  undetermined)  functional  of  the  values 
of  all  the  Fourier  amplitudes.     If  now  ^  represents  the  direct  interaction  terms 
and  appeal  is  made  to  the  weak  dependence  hypothesis,   the  ensemble  averages  on 
the  right  side  of  (2,U)  may  be  reduced  to  expressions  involving  only  the  ensemble 
averages  ^S^.  .(kjt,s)^,   ...  of  the  impulse  response  tensor  and  covariances  of 
the  form  <,u,  (k,t)u     (k,s)^,    .    •   .      .To  complete  the  evaluation  of  the  triple 
moments  it  remains  to  determine  the  averaged  impulse  response  tensor.     Equations 
of  motion  for  this  tensor  may  be  obtained  by  considering  further  the  response  of 
the  system  to  arbitrary  external  perturbations.     These  equations  and  the  equations 
of  motion  for  the  covariance  tensor  serve  to  determine  the  tira  tensors  simul- 
taneously.    The  complete  derivation  for  the  stationary  case  is  given  in  Paper  I. 

Using  an  extension  of  the  direct  interaction  hypothesis,   certain  fourth- 
order  moments  of  the  Fourier  amplitudes  also  may  be  expressed  in  terms  of  the 
covariance  tensor  and  the  averaged  impulse  response  tensor,  and  a  corresponding 
reduction  may  be  made  for  certain  third-order  moments  involving  random  external 
forces  which  may  be  driving  the  system.     It  is  shown  in  Paper  I  that  from 
knowledge  of  these  moments  it  is  possible  to  obtain  the  complete  set  of  higher- 
order  moments  of  a  distribution  which  satisfies  the  equations  of  motion  exactly, 
aiid  thus  to  obtain  a  complete  statistical  description  of  the  turbulence.     In  the 
present  paper,  however,  we  shall  be  concerned  only  with  second-  and  third-order 
moments. 
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3.     FORMUL/STION  PDR  STATIONARY  ISOTPDFTC   TUEBULENCE 

The  theorj'  described  above  takes  its  simplest  form  for  the  stationary 
isotropic  case.     The  concept  of  stationary  isotropic  turbulence  involves  a 
considerable  idealization  of  physical  flows.     For  turbvilence  to  be  stationary 
there  must  be  driving  forces  to  replace  energy  lost  by  viscosity;   in  ortler  to 
maintain  isotropy  the  driving  forces  must  themselves  be  isotropic,  and  this  is 
physically  unrealistic.      (We  may  visualize  such  forces,  perhaps,   as  due  to  a 
random  volume  distribution  of  stirring  devices.)     Despite  the  unrealistic  nature 
of  stationary  isotropic  turbulence,  it  is  valuable  for  study.     Analytically,   it 
is   the  simplest  case;  physically,   it  is  plausible,   if  the  driving  forces  invoked 
are  confined  to  low  wavenumbers,   that  the  structure  at  high  wavenumbers  will  not 
differ  appreciably  from  that  of  freely  decaying  isotropic  turbulence   (also  an 
idealization).     We  shall  return  to  this  question  after  exploring  the  consequences 
of  the  stationary  isotropic  theory. 

The  stationary  isotropic  state  may  be  described  in  a  mathematically 
consistent  fashion  by  replacing  (2.2)  with 

(3.1)  u.(k,t)  +  vk^u.(k,t)  -  -  ilc  P..(k)     Z      u  (p,t)u  (q,t)  +  f.(k,t), 

where  f.(k,t)   is  a  statistically  stationary  and  isotropic  solenoidal  forcing 
term.     An  alternate  procedure,  which  in  some  ways  is  more  closely  related  to  the 
physical  mechanism  of  support  by  a  mean  shear  flow,  is   the  introduction  of 
negative  viscosity  at  low  wavenumbers  rather  than  external  forcing  terms. 

The  conditions  of  stationarity  and  isotropy  require  that  the  covariances 
constructed  from  u(k,t)  and  f(k,t)  have  the  forms 
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(3.2)  (L/2n)\u^(k,t  +r)uj*(k«,t)>   =   ^  5_,j^P^^(k)U(k,T), 
(L/2n)\f^(k,t  +T)fj*(k',t)>   =   I  5,^P^j(k)F(k,r), 

(L/2n)^Re<f^(k,t  +  r  )u^*(k',t)>  -   |  5,^,P^^(k)G(k,T), 

where  the  scalars  U,  F,  and  G  are  all  real  and  U  and  F  are  even  functions  of  -c 
In  the  limit  L->  oo    (which  is  required  for  rigorous  isotropy), 

(3.3)  U(k,r)  -  (2tt)"^  I  U(x,T)e"^-'^d\, 

where  U(x,t)  = -(u.  (x  +  y,  t  +t)u.  (y,t)),  with  corresponding  expressions  for 
F(k,T)  and  G(k,T).* 

From  (3.1)  we  may  obtain  the  equation  of  motion 

(3.b)         U(k,-c)  +  vk^U(k,r)  =  S(k,T)  +  G(k,T), 
where  the  dot  now  signifies  -r  differentiation  and 

(3.5)         S(k,T)  -  (L/2n)^  Im{k  X  <u  (p,t)u  (q,t)u  *(k,t- r  )>]  . 

The  perturbation  procedure  described  in  the  preceding  section  may  be  used  to 

reduce  the  triple  moment  on  the  right  side  of  (3.5).  Since  the  turbulence  is 

stationary  and  isotropic,  the  ensemble  average  of  the  impulse  response  tensor 

must  display  these  properties  also.  Thus,  (2,5)  leads  to 

In  accordance  with  the  convention  adopted  in  Paper  I,  the  normalization  in 
(3.2)  corresponds  to  the  requirement  that  the  velocity  obey  cyclic  boundary 
conditions  on  opposite  faces  of  a  cube  of  side  L,  not  necessarily  vanish. 
To  avoid  confusion  with  another  notation,  f . (k,t)  was  denoted  by  F^(k,t)  in 
Paper  I.  F(k,T)  has  the  same  meaning  in  both  papers. 
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(3.6)  6u^(k,t)  -  j     g(k,t-3)C^(k,3)d3, 

where  6  denotes  the  ensemble  averaged  induced  perturbation  and  the  Impulse 
response  function  g(k,t-3)   is  defined  by  P^,(k)g(k,t-s)  = <i^.  .(kjtjs)).     Then, 
as  detailed  in  Paper  I.   the  perturbation  procedure  yields,   in  the  limit  L-*-  °^  , 

(3.7)  S(k,T)  =  nkjj    pqdpdqj^    d3[a(k,p,q)g(k,s)U(p,3  +T)u(q,s  +  t) 

-b(k,p,q)g(p,s)n(k,s  --r)U(q,s)]    , 

where  the  domain  of  integration  extends  over  all  wavenumbers  p  and  q  such  that 
p,q,  and  k  can  form  the  legs  of  a  triangle.  The  quantities  a  and  b  are  defined 

by 

(3.8)  a(k,p,q)  =  (l-xyz-2y^z^). 


—1        2   2    3     2 
b(k,p,q)  =■  2k  p(-z+xy+x  z+y  z+2z  +2x7z  ), 


where  x,  y,  z  are  the  cosines  of  the  interior  angles  opposite  the  legs  k,  p,  q, 

respectively.  They  obey  the  identities 

(3.9)  a(k,p,q)  «  a(k,q,p)  >  0, 

k  b(k,p,q)  «  p'T)(p,k,q), 

b(k,p,q)+b(k,q,p)  =«  2a(k,p,q), 

which,  as  we  shall  see  shortly,  are  connected  with  the  energy  conservation 
properties  of  the  interaction. 

Since  the  driving  forces  are  arbitrary  in  any  event,  let  us  take  them 
to  be  Gaussianly  distributed.  Then  they  may  be  regarded  as  sums  of  infinite 
numbers  of  infinitesimal  statistically  independent  increments.  Each  increment 
individually  may  be  considered  an  infinitesimal  perturbation  to  the  equations  of 
motion  J  we  thereby  find 
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(3.10)  G(k,T)  =  J  g(k,3)F(k,s  +-r)ds. 

By  considering  the  effect  of  an  arbitrary,  statistically  stationary 
and  isotropic,  solenoidal  perturbation  6f. (k,t)  on  a  statistical  equation 
obtained  by  multiplying  (3»l)  through  by  f.  (k,t')  and  averaging,  we  may  obtain 
the  equation  of  motion 

(3.11)  g(k,x)  +  vk  g(k,T)  =  -Ttkjj    pqdpdqb(k,p,q) 


X  )  g(k,T-s)g(p,3)U(q,3)ds     {-^>  0), 


as  described  in  Paper  I,  For  t<  0,  g(k,T)  =  0,  since  there  can  be  no  response 
to  a  perturbation  before  it  is  applied.  If  F(k,-r)  is  prescribed,  (3.U),  (3.7) » 
(3»10),  and  (3.11)  form  a  complete  system  determining  U(k,T)  and  g(k,T),  provided 
suitable  boundary  conditions  are  imposed.  The  condition  of  static narity  requires 

(3.12)  U(k,0)  =  0, 

which  may  be  formulated  by  (3.U)  and  (3.7)  as  an  integral  condition  (U.2) 
expressing  the  balance  of  energy.  The  other  required  boundary  condition  is 

(3.13)  g(k,  +  0)  -  1 

which  follows  immediately  from  the  definition  of  g(k,T)  and  the  fact  that  the 
equation  of  motion  (3.1)  is  first-order* 

U.     THE  SPECTRAL  TR/.NSPORT  OF  ENERGY 

The  covariance  scalars  introduced  in  the  last  Section  may  be  written 

(U.l)  \  U(k,T)  -  (Unk^)"^(k)r(k,T), 

F(k,T)  -  (UTtk^)"-'T(k)ii(k,T), 
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where  E(k)   is  the  kinetic  energy  spectrum  and  F(k)  the  forcing  spectrum,  both 
per  unit  mass  of  fluid.     The  total  kinetic  energy  per  unit  mass  is   j       E(k)dk, 
The  quantity  r(k,T)  is  the  time  conflation  function  of  the  Fourier  mode, 
normalized  so  that  r(k,0)  =  1,  with  a  corresponding  interpretation  for  n(k,'c). 
If  we  rewrite   (3«U)   in  terras  of  the  new  quantities,  set  t  =  0,  and  note  the 
symmetry  of  the  integration  domain  with  respect  to  p  and  q,  we  obtain  the 
relation 

(U.2)  2vk^E(k)  =  III  S(k|p,q)dpdq  +  F(k)£   g(k,3)ti(k,s)d3, 

where 

(U.3)  S(k|p,q)  =  (k/pq){2k^a(k,p,q)e(k,p,q)E(p)E(q) 

-[p\(k,p,q)e(p,k,q)E(q)  ♦  q^b(k,q,p)  ©  (q,k,p)E(p)]E(k)} 

and 

(U.U)  e(k,p,q)  «   j     g(k,s)r(p,s)r(q,3)ds, 

■'o 

The  left  side  of  (U«2)  is  the  mean  power  dissipated  by  viscosity,  per 
unit  mass  and  per  unit  wavenumber.  The  second  term  on  the  right  may  be  inter- 
preted as  the  mean  power  input  by  the  driving  forces,  per  unit  mass  and  per  unit 
wavenumber.  (Hereafter,  all  powers  will  be  understood  to  refer  to  unit  mass.) 
It  should  be  noted  that  the  linear  dependence  of  this  term  on  F(k)  is  only 
apparent  J  the  response  function  g(k,8)  is  an  implicit  function  of   and  \^ 
determined  by  (3.U),  (3. ID),  and  (3,11),  The  first  term  on  the  right  may  be 
interpreted  as  a  mean  power  input  due  to  nonlinear  interaction  of  the  velocity 
modes.  Thus,  S(k|p,q)dkdpdq  represents  the  mean  power  delivered  to  the  Fourier 
modes  in  the  interval  dk  as  a  consequence  of  their  triad  interactions  with  all 
pairs  of  modes  of  which  one  member  lies  in  each  of  the  intervals  dp  and  dq. 
S(k|p,q)  is  syTTTmetric  in  p  and  q,  and  the  factor  w  in  (li.2)  occurs  because  each 
pair  of  modes  is  counted  twice  in  the  integration. 
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lt may  be  deduced  readily  from  (3.1)  that  the  interaction  of  any 
tried  of  Fourier  modes  is  individually  conservative j  this  implies 

(U.?)  S(klp,q)  +  S(p|q,k)  +  S(q|k,p)  =  0 

which  may  also  be  obtained  easily  from  (U.3)  by  using  the  trigonometric 
identities  (3.9)  and  the  fact  that  9(k,p,q)  is  sjonmetric  in  p  and  q.  The  overall 
conservation  property 

(U.6)  i  f  dk  I f  S(klp,q)dpdq  -   0 

^  Jq        j) c^ 

follows  directly  from  the  detailed  conservation  relation. 

It  should  be  noted  that  while  the  concept  of  the  elementary  interaction 
of  a  triad  of  Fourier  modes  has  a  particular  significance  because  of  the  detailed 
conservation  property,  the  same  is  not  true  of  the  concept  of  exchange  of  energy 
between  pairs  of  modes,  which  has  been  used  frequently  in  turbulence  theory. 
There  is  no  way  to  tell  how  much  of  the  energy  a  mode  k  receives  from  interaction 
with  modes  p  and  q  comes  from  p  rather  than  from  q  (Paper  I,  footnote  7). 

The  quantities  which  enter  expression  (U.3)  may  be  given  fairly  direct 
physical  interpretations.  We  shall  consider  first  the  factors  9(k,p,q),  e(p,k,q), 
9(q,k,p)  which,  from  their  dimensions,  evidently  represent  characteristic  times 
for  the  triad  interaction.  We  may  anticipate  that  these  factors  ai^  essentially 
positive.  The  physical  significance  of  the  response  function  g(k,s)  is  that  it 
traces  the  relaxation  of  an  initial  applied  excitation  of  mode  k  through  the 
joint  action  of  viscous  decay  and  the  energy  exchange,  or  mixing  action,  with 
the  infinitely  large  number  of  only  weakly  statistically  inter-related  other 
modes.  We  may  expect  that  it  falls  to  zero  with  increasing  s  in  a  fairly  smooth 
fashion  and  that  it  does  not  exhibit  strongly  negative  regions.  The  correlation 
function  r(k,s)  also  is  determined  by  the  general  mixing  action.  We  may  expect 
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that  it,  too,  is  essentially  positive  and  falls  fairly  smoothly  to  zero  with 
increase  of  s.  Thus  ©(k,p,q),  ©(p,k,q),  and  9(q,k,p)  may  be  expected  to  be  the 
order  of  characteristic  correlation  and  relaxation  times  of  the  three  modes 
involved. 

The  proportionality  of  the  terms  in  S(k|p,q)  to  the  6  factors  may  be 
interpreted  as  follows.  A  non-zero  mean  transfer  of  energy  among  the  three 
modes  inquires  that  certain  phase  relations  among  these  modes  be  established 
in  a  statistical  sense.  These  relations  are  built  up  by  the  direct  interaction 
of  the  three  modes,  according  to  the  direct  interaction  hypothesis,  but  are 
simultaneously  broken  down  by  the  viscous  decay  and  effectively  random  mixing 
with  all  the  other  modes,  as  expressed  by  the  g  and  r  functions.  Loosely 
speaking,  we  may  say  that  9(k,p,q)  is  an  effective  time  during  which  the  triple 
phase  correlation  can  build  up  before  it  is  broken  down. 

The  first  of  identities  (3»?)  indicates  that  a(k,p,q)  is  never  negative, 
and  the  third  that  b(k,p,q)  is  typically  positive.  If,  as  anticipated,  the  © 
factors  are  positive,  the  terms  in  {k»3)   involving  a  represents  a  positive  flow 
of  energy  to  mode  k  while  those  involving  b  represent  a  typically  negative  flow. 
The  net  flow  is  the  resultant  of  these  absorption  and  emission  terms.  It  will 
be  noticed  that  in  contrast  to  the  absorption  term  the  emission  terms  are 
proportional  to  E(k),  This  indicates  that  the  energy  exchange  acts  to  maintain 
equilibrium.  If  the  spectrum  level  were  suddenly  raised  to  much  higher  than 
the  equilibrium  value  in  a  narrow  neighborhood  of  k,  the  emission  terms  would 
be  greatly  increased  while  the  absorption  terra  would  be  little  affected^  thus, 
energy  would  be  drained  from  the  neighborhood  and  equilibrium  re-established. 
The  structure  of  the  emission  and  absorption  terms  is  such  that  in  general  we 
may  expect  the  eneri^   flow  to  be  from  strongly  to  weakly  excited  modes,  in 
accord  with  general  statistical  mechanical  principles. 
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Although,  as  was  noted  above,  it  is  impossible  to  define  unambiguously 
the  transfer  of  energy  between  a  pair  of  modes,  it  is  meaningful  to  define  a 
transport  power  Il(k)  as  the  mean  power  input  to  all  modes  of  waveramber  higher 
than  some  value  k  from  all  modes  of  wavenumber  less  than  k.     For  this  purpose 
we  divide  all  the  triad  interactions  into  four  classes  as  shown  in  Figure  1, 
Classes  A  and  B  do  not  transfer  energy  across  the  boundary.     II(k)  is  therefore 
the  power  input  to  all  modes  k'  >  k  by  all  interactions  in  class  C  less  the 
power  input  to  all  modes  k'<  k  by  all  interactions  in  class  D,     Thus, 

(U.7)  TKk)  =1/    dk'IJ  ,     S(k'|p,q)dpdq  -  |/   dk'||^,      S(k' |p,q)dpdq. 

The  upper  limit  on  the  first  integral  sign  is  taken  as  2k  instead  of  oo  because 
for  k'>  2k  there  are  no  p  and  q  such  that  p,q  <  k  and  k',  p,  q  form  a  triangle. 

It  follows  from  the  conservation  properties  that  the  interaction  power- 
input  density  appearing  in  (U.2)  is  related  to  IT(k)  by 

(U.8)  i  //  S(klp,q)dpdq  -  -  ^^ 


5.     RESIDNSE  AND  TIME  CORRELATION  FUNCnONS  FOR  HIGH  WAVEKTJMBERS 

5»1»     The  response  function 

The  equations  of  motion  (3.1j)  and  (3.11),  with  S(k,T)  and  G(k,T)  given 
l>y  (3*7)  and  (3,3D),   are  of  a  type  which  does  not  seem  to  have  been  studied,  and 
it  may  be  anticipated  that  in  general  their  solution  will  present  severe  diffi- 
culties.    It  is  not  too  difficult,  however,   to  develop  a  solution  which  is 
asymptotically  valid  for  the  inertial  and  dissipation  wavenumber  ranges  in  the 
limit  of  very  high  Reynolds  number.     The  simplest  starting  point  is  the  equation 
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of  motion  for  g(k,T). 

Using  (U.l)  we  may  rewrite  (3.11)  in  the  form 

(5.1)   g(k,T)  +  vk^g(k,T)  >-  |J/dpdq£b(k,p,q)E(q)r  g(k,T-s)g(p,s)r(q,s)ds 

(t:>0). 

Now  let  us  consider  a  k  sufficiently  high  that  all  but  a  very  small  fraction  of 
the  total  energy  lies  very  far  below  k.  Because  of  the  factor  E(q),  it  appears 
plausible  to  assume  that  the  right  side  of  (5.1)  should  be  dominated  by  contri- 
butions from  the  principal  energy  containing  region  q  «  k.  This  assumption  is 
supported  by  direct  inspection  of  the  original  equation  of  motion  (3.1).  It  will 

be  noticed  that  the  coefficient  ik  P. .(k)  of  the  bilinear  terms  is  independent 

m  ij  - 

of  IpI  and  jq].  Apart  from  certain  angular  dependence  effects,  the  coupling  of 

the  mode  k  to  modes  p  or  q  «  k  appears  quite  as  strong  intrinsically  as  its 

coupling  to  modes  -^k.  Thus  we  might  expect  that  the  contribution  of  the 

bilinear  terms  to  the  motion  should  be  dominated  by  terms  involving  the  modes 

which  contain  most  of  the  total  excitation.  Also,  it  is  clear  from  the  appearance 

of  the  factor  k  that,  in  general,  the  rate  of  change  of  u.(k,t)  ascribable  to 
m  1 

the  nonlinear  terms  tends  to  increase  with  k.  Thus  we  may  anticipate  that  the 
characteristic  times  of  the  functions  g(k,s)  and  r(k,s)  decrease  with  increasing 
k. 

Our  assumptions  are  related  to  the  observations  in  the  x  space  representa- 
tion that  the  time  variation  of  the  fine  structure,  viewed  in  a  fijced  coordinate 
system,  should  be  dominated  by  the  sweeping  action  of  the  energetic  large  scale 
motion  and  that  small  structures  should  be  swept  past  fixed  points  more  rapidly 
than  large  ones. 


-  13  - 


Returning  to  (5.1),  we  note  that  for  q  «  k  the  triangle  relation  gives 
p  =s  kj  then  we  may  approximate  g(p,s)  by  g(k,3),  provided  it  is  a  reasonably 
smooth  function.  Since  the  time  characterizing  r(q,s)  should  be  very  large 
compared  to  that  of  g(k,s),  we  shall  assume  that  in  the  region  where  g(k,3)  is 
significantly  different  from  zero  we  may  replace  r(q,3)  by  r(q,0)  =  1.  Finally, 
we  may  find  from  the  definition  (3.8)  that,  for  k  ~  p, 

2 
b(k,p,q)  ~  2sin  p, 

where  p  is  the  interior  angle  between  k  and  q.  The  wavenumber  integration  for 
q  in  the  energy  range  now  can  be  readily  carried  out,  yielding 

■1 


pq"-H)(k,p,q)E(q)dpdq  *  2   dqE(q)J    pq'-'-sin^pdp-'?  k   E(q)dq, 


where  k  «  k  is  a  wavenumber  below  which  lies  very  nearly  all  the  energy.  Then, 

noting 

■  e  r°° 

E(q)dq    ~  E(q)dq  =4v   ^, 

where  v     is  the  rras  velocity  in  any  direction,  we  obtain  the  equation  of  motion 

(5.2)  g(k,T)  +  vk^g(k,T)  =  -v^VJ     g(k,T-s)g(k,3)ds.  (t>  0). 

If  our  assiunptions  and  approximations  are  well-founded,   this  equation  should  be 
asymptotically  exact  for  sufficiently  high  k. 

The  solution  of  the  Laplace  transform  of  (5.2)  for  the  boundary  con- 
dition (3.13)  is  a  standard  form,  and  we  thereby  find 

(5.3)  g(k,T)  -  exp(-vk^T)J^(2v^kT)/(T^kT)  (t>0). 
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A  detailed  investigation  of  the  behavior  of  the  right  side  of  (5.1) >  with  g 
given  by  (5.3)  and  r  estimated  from  the  results  to  be  derived  in  Section  5.2, 
indicates  the  consistency  of  our  assumptions  and  approximations.  It  is  found 

that  the  restriction  of  the  integration  region  to  q  «  k  can  give  rise  to 

2  —1 

appreciable  fractional  error  in  g(k,T)  only  when  ^  »  (vk  )  .  For  such  i:  , 

error  can  arise  because  the  integral  over  s  may  be  much  larger  for  p,  q  --  ik  than 
for  p  »  k.  However,  g(k,T)  is  very  small  for  such  x  ,  and  the  possible  error 
is  unimportant.  Moreover,  it  decreases  with  increasing  separation  of  k  and  the 
energy  containing  region.  Our  further  approximations  are  found  to  lead  to  no 

difficulties. 

2 
The  factor  exp(-vk  t)  in  (5.3)  is  the  response  function  appropriate  to 

the  decay  of  an  applied  excitation  by  viscosity  without  interaction  with  the 

other  modes.  The  factor  J,(2v  kT)/(v  kx)  (plotted  in  Figure  2)  represents  the 

relaxation  associated  with  the  energy  exchange  between  the  given  mode  and  others 

with  very  nearly  equal  wavevectors,  due  to  the  mixing  action  of  the  energy 

containing  region.  By  this  process,  the  applied  excitation  is  mixed,  or 

diffused,  into  many  neighboring  (but  statistically  weakly  dependent)  modes. 

The  characteristic  time  (v  k)"  associated  with  the  mixing  contribution  to  g(k,T) 

is  the  same  order  as  that  which  would  be  expected  from  the  convection  of  a 

periodic  pattern  of  wavenumber  k  by  a  uniform  velocity  of  magnitude  v  .  However, 

the  action  of  the  low  wavenuitibers  in  the  present  theory  is  very  different  from 

the  convective  action  of  a  uniform  field.  Such  a  field  does  not  produce  mixing 

of  neighboring  modes  and  therefore  gives  rise  to  entirely  coherent  time  dependence, 

rather  than  relaxation  effects.  This  matter  will  be  discussed  further  in 

Section  9.1. 
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In  accord  with  earlier  anticipations,   the  response  function  given  by 
(5»3)  falls  off  fairly  smoothly  with  increase  of    r.     However,   there  are 
oscillations,   as  shown  in  Figure  2.     Although  these  are  small,   and  therefore  do 
not  have  much  affect  on  the  9  factors  and  the  magnitude  of  energy  transport, 
they  have  a  considerable  intrinsic  interest.     In  the  language  of  linear  systems, 
they  indicate  that  the  mixing  with  neighboring  modes  presents  a  partly  reactive 
impedance  to  a  given  mode.     This  may  be  understood  further  in  connection  with  a 
physical  interpretation  of  the  structure  of  the  right  side  of  (5.1).     The 
transfer  of  part  of  an  initially  applied  excitation  from  mode  k  to  mode  p 
involves  two  steps.     First,  a  perturbation  is  induced  in  mode  p.     This  is 
determined  by  the  response  function  g(p,3),  and  involves  also   the  amplitude  of 
a  third  mode,   q,  in  accord  with  the  nonlinear  structure  of  (3.1).     Second,   the 
perturbation  in  mode  p  reacts  upon  mode  k  producing  in  it  a  change  determined  by 
the  response  function  g(k,T-s),  and  again  involving  the  amplitude  of  mode  q. 
The  role  of  mode  q  in  the  two  steps  accounts  for  the  expression  E(q)r(q,s).     It 
is  clear  that     for  part  of  the  energy  of  the  initial  excitation  to  be  removed 
from  mode  k     the  reaction-induced  perturbation  must  be  effectively  of  opposite 
phase  from  the  initial  disturbance.     The  oscillations  in  g  indicate  that  at 
certain  difference- times  the  counter-perturbation  actually   "overshoots'  so  that 
there  is  in  effect  a  small  coherent  oscillation  of  the  initial  excitation  energy 
between  modes  k  and  p. 

It  should  be  noted  that  g(k,T)  plays  a  double  role.     It  gives  the 
(average)  response  to  individual  perturbations,  but  in  doing  so  it  must  take 
account  of  the  dynamical  effect  of  an  infinite  number  of  repetitions  of  indivi- 
dual pairs  of  actions  and  reactions  as  just  described.     This  is  reflected  in  the 
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inhomogeneity  of  (5.1)   and  (5.2)  in  g.     The  right  sides  may  be  expanded  by 
iteration  into  forms  involving  all  powers  of  g,   and  many-fold  tinve  integrals. 
This  reason  for  the  nonlinearity  of  (5,1)   (which  is  of  higher  degree  than  that 
of  (3.1))  would  exist  even  if  the  equations  of  motion  of  the  system  were  linear. 

5.2«     The  time  correlation  function 

Measurements  on  real  turbulent  flows  seem  to  indicate  that  the  supply  of 
energy  from  the  mean  flow,  or  from  whatever  other  external  source  may  be  acting, 
is  effectively  confined  to  wavenumbers  in  and  near  the  energy  containing  region. 
Then  we  may  rewrite  (3.1l),  with  S(k,T)  given  by  (3.7),   in  the  form 

(5.U)  r(k,T)  +  vk2r(k,T)  -|J        k^a(k,p,q)^^PJffiM    g(k,3-T)r(p,3)r(q,s)d3 


A 


-p^(k,p,q)E(q)/     g(p,3)r(k,T-s)r(q,s)d3l  ^  . 


(In  obtaining  this  form  a  slight  change  of  time  variable  has  been  made,  and  it 

has  been  noted  that  r  is  an  even  function  of  time.) 

Now  let  us  assiune  that  the  Rejmolds  number  is  high  enough  that  there  is 

a  range  of  k  large  compared  to  energy  containing  wavenumbers  but  small  compared 

to  T  /v,  so  that  vk  «  v  k.  For  such  k  we  may  anticipate  that  the  behavior  of 
0  o 

r(k,T)  should  be  dominated  by  the  energy  range  mixing  action  discussed  previously 
and  consequently  be  characterized  by  a  time  the  order  of  (v  k)  .  If  this  be  so, 
the  term  vk  r(k,T)  on  the  left  of  (5.U)  should  be  negligible  compared  to  r(k,T), 
except  in  the  immediate  vicinity  of  the  origin.  Now,  following  arguments  similar 
to  those  used  in  obtaining  (5.2),  let  us  retain  on  the  right  of  (5.U)  only  those 
contributions  involving  b  such  that  E(q)  refers  to  the  energy  containing  range 
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and  only  those  contributions  involving  a  such  that  E(p)  or  E(q)  refer  to  this 
range.  Making  the  further  approximations  used  in  obtaining  (5.2),  taking  either 
E(q)/fe(k)  or  E(p)/^(k)  ^    1,  as  appropriate,  using  the  fact  that 


and  writing 


2 

a(k,p,q)  ~    sin  p      (q  «  k). 


ds  "  (   ds'  +  I   ds,  where  s'  =x-s, 
r 


we  obtain  the  asymptotic  equation 

(5.5)     r(k,T)  -  -T^VJ  [   g(k,T-s')r(k,s')ds' 


[r(k,T-s)g(k,3)  -  g(k,s-T)r(k,s)]  dsV  . 

The  solution  satisfying  (3.12)  is  r(k,-c)  «  g(k,  |x|),  since  for  this 

choice  the  second  integral  on  the  right  of  (5.5)  vanishes  and  the  equation 

2 
reduces  to  the  asymptotic  form  taken  by  (5.2)  for  the  case  vk  «  v^k.  Thus 

we  have  the  result 

(5.6)  r(k,T)  -  J^(2v^k-c)/(v^kT). 

A'  consistency  check  similar  to  that  described  after  (5.3)  shows  that  this 
solution  is  indeed  asymptotically  valid  if  k  satisfies  the  double  inequality 
we  have  postulated.  The  frequency  spectrum  corresponding  to  our  solution  is 

_r     /    \2t1/2 

(5.7)  9(k,co)=|[l.  (^)J        (co<2y^k), 

-  0  (oo  >  2v  k). 

This  function  is  shown  in  Figure  3. 
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The  spectrum  of  a  correlation  function  must  be  positive-definite,  and 
the  fact  that  the  r(k,a3)  given  above  has  this  property  is  some  evidence  of  the 
consistency  of  our  theory  and  of  the  realizability  of  solution-distributions 
obeying  our  statistical  hypotheses.  The  sharp  cut-off  in  the  spectrum  at 
CO  =  2v  k  is  a  consequence  of  the  asymptotic  nature  of  our  results.  We  have 
explicitly  neglected  the  time  variation  of  mode  k  due  to  interactions  with  modes 
of  wavenumber  appreciably  greater  than  k.  At  any  high  but  finite  k,  the  inter- 
action with  higher  wavenumbers,  having  higher  characteristic  frequencies,  must 
contribute  a  slight  tail  to  the  spectrum,  extending  to  infinite  frequency.  The 
height  of  the  tail  decreases  with  increasing  separation  of  k  from  the  energy 
containing  region, 

6.  THE  INERTIAl  RAMiE 

Let  us  assume  the  existence  of  an  extended  range  of  k  which  satisfies 
all  the  inequalities  invoked  in  Section  5  and  for  which  the  dissipation  term  on 
the  left  of  (U»2)  is  negligible.  The  energy  balance  equation  then  is  simply 


(6.1)  |||^S(klp,q)dpdq  -  0. 


It  follows  according  to  the  analysis  of  Section  5  that  the  contributions  to  this 
equation  for  which  p  or  q  is  in  the  energy  containing  range  are  given  asymptoti- 
cally by  2E(k)  times  the  right  side  of  (5.5)  for  x  »  0.  The  s  integrals  in  this 
expression  vanish  identically  (independently  of  the  form  of  r(k,s)  and  g(k,s)), 
which  suggests  that  modes  of  high  wavenumbers  do  not  receive  energy  from  direct 
interaction  with  the  energy  containing  range.  This  is  in  accord  with  the  widely 
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accepted  intuitive  notion  that  the  transfer  of  energy  up  the  spectrum  proceeds 
by  an  essentially  local  cascade  process.  We  shall  assume  that  this  is  so, 
determine  the  spectrum  in  the  inertial  range,  and  then  check  the  consistency  of 
our  assumption  by  obtaining  an  expression  for  the  energy  transport  due  to  inter- 
action among  distant  wavenumbers. 

The  existence  of  an  effectively  local  cascade  process,  together  with 
the  conservative  character  of  the  interaction,  implies  that  in  the  inertial  range, 
where  dissipation  is  negligible,  the  transport  power  must  satisfy 

(6.2)  n(k)  =  e 

where  the  parameter  e  is  the  mean  power  eventually  dissipated  at  very  high 
waveniimbers.  In  accord  v;ith  the  concept  of  a  local  cascade  process  we  shall 
assvime  that  contributions  to  IKk),  as  given  by  (U.7)  and  (lio3)>  may  be  neglected 
if  p,  q,  or  k'  is  verj"^  large  or  small  compared  to  k.  For  all  the  remaining 
contributions,  the  relevant  response  functions  and  correlation  functions  are  of 
the  form  J,(2v  kT)/(v  kx).  The  9  factors  defined  by  (U.li)  then  are  symmetric  in 
their  arguments  and  have  the  dimensional  form 

[e(k,p,q)]  -[vj-^[k]-^ 
Consequently  XT(k)  has  the  dimensional  fonn 

[n(k)]  =  [vi-^[k]3[E(k)j2 


It  seems  rather  clear  from  this  that  if  n(k)  is  to  be  determined  wholly  by 
contributions  with  k',  p,  and  q  all  in  some  essentially  local  neighborhood  of  k 
and  is  to  have  the  value  e,  which  is  independent  of  k,  then  one  must  have 
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(6.3)  B(k)  -  f(0)(e^r^)^/V^/^ 

where  f(0)  is  a  numerical  constant.  (We  choose  this  symbol  for  later  convenience.) 

The  argument  just  given  only  suggests  the  necessary  form  of  E(k),  and 
we  must  check  to  see  that  (6.3)  actually  satisfies  (6.1)  and  leads  to  a  local 
cascade  process.  If  (6.3)  is  substituted  in  (U.3)  and  one  takes  g(k,s)  =  g(v  ks), 
r(k,s)  «  r(v  ks)  (without  specifying  the  particular  functional  forms  of  r  and  g), 
it  is  not  difficult  to  verify  by  formal  manipulations,  using  the  identities  (3«9), 
that  (6.1)  is  satisfied.  This  actually  is  rather  clearly  implied  by  the  conserva- 
tive nature  of  the  theory  and  the  dimensional  considerations  above.  The  formal 
property  has  meaning,  however,  only  if  the  integrals  involved  converge  properly, 
which  corresponds  physically  to  the  presumed  localness  of  cascade.  We  may  verify 
the  latter  property  by  considering  the  total  power  input  to  all  modes  of  wave- 
number  k"  >  k'  f3X»m  direct  interactions  with  all  mode  pairs  p,  q  such  that 
p  or  q  <  k  «  k',  where  k  and  k'  are  fixed  wavenunibers  and  all  waverairabers  con- 
cerned are  within  the  inertial  range.  Using  r(k,s)  ■=  g(k, |s|)  "  J^(2v^ks)/(v^ks) 
and  computing  this  power  as 

(por  q<K) 
.00 


(6.1;)  n(k'|k)=J   dk"       S(k"|p,q)dpdq, 

we  find,  asymptotically, 

(6.5)         n(k'|k)  =  (numerical  factor)  [f (O)l^e(kA')^'^     (k  «  k'). 

The  triad  interactions  involved  in  the  integration  are  shown  in  Figure  U. 
Tl(k'|k)  goes  to  zero  with  kA'>  and  the  cascade  process  is  asymptotically  local, 
9s  originally  assumed.  The  dependence  on  kA'  is  not  particularly  strong, 
however,  and  the  cascade  is  rather  ill-defined. 
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It  remains  to  determine  the  c»nstant  f(0).     This  may  be  done  by 
substituting  (6,3)  into  (1;.7)»  taking  the  form  J,(2r  ks)/(v  ks)  for  the  response 
and  correlation  functionsj  evaluating  the  definite  multiple  integrals,   and 
noting  (6,2).     After  this  is  done,   (5.6)  and  (6,3)  provide  a  complete  asymptotic 
solution  for  the  inertial  range. 

If  the  inertial  range  extends  through  and  well  above  the  wavenuniber  k, 
the  total  energy  above  k  according  to   (6.3)  is 


(6.6)  f   E(p)dp  «  if  (0)(e7^)^/V^/l 


Now  let  us  consider  the  extrapolation  of  this  expression  down  to  low  k.  Assuming 

that  f(0)  is  crudely  the  order  of  unity,  it  is  clear  that  (6,6)  gives  an  energy 

2 
the  order  of  v   (per  unit  mass)  when  k  is  the  order  of 


(6.7)  k^  =  eA^^. 

Such  an  extrapolation  violates  our  conditions  for  the  inertial  range,  of  course, 
but  nevertheless  this  result  suggests  strongly  that  k  is  actually  a  wavenumber 
characteristic  of  the  energy  containing  range.  It  is  also  known  empirically 
that  this  is  the  case(Batchelor,  1953,  p.l03).  We  therefore  may  define  a 
Reynolds  number  characteristic  of  the  energy  containing  range  by 

(6.8)  R^  =  v^k^-Vv  =  v^Vev. 

If  the  inertial  range  extends  through  and  well  below  a  wavenumber  k, 
the  total  power  dissipated  below  k  is 
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(6.9)  j    2vp2E(p)dp  -  ^f(0)v(eT^)^/V/2. 

This  expression  becomes  the  order  of  e  for  k  the  order  of 


(6.10)  k,  .  e^^^v  -^/V2/3  -  R  2/3k  . 

do  o        o 


It  ia  clear,  then,  that  necessary  conditions  for  the  validity  of  our  asymptotic 
inertial  range  spectntm  are 


(6.11)  k  «  k  «  YL^'\^, 

for  an  extended  range  of  k.  This  requires  in  turn  that  R  satisfy 

(6.12)  R^^/3  »  1. 

Since  v  /v  »  R  k  ,  it  is  clear  that  the  condition  vk  «  ▼  k,  invoked  in  the 
o     o  o'  o  ' 

derivation  of  (5.6),  is  contained  in  (6.11). 

A  Reynolds  nuinber  R,  may  be  defined  by  R.  «  v  k- "  /w,  where 
t^  -  eASv^^v  (Batchelor,  1953,  p.^D.  Thus, 


(6.13)  \   -  (15RQ)^^^ 


The  condition  (6.12)  is  equivalent  to 
(6.II1)  R^^/3  »  1. 
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7.  THE  DISSIPAHON  RANGE  A.T  HUGE  REYNOLDS  NUMBERS 


7,1,  The  integral  equation  for  E(k) 


The  discussion  at  the  end  of  Section  6  indicated  that  a  major  part  of 

3  p] 
1/3 


2/3 

the  dissipation  takes  place  at  wavenuiribers  the  order  of  k,  =  R        k  ,     This 

a.    0    0 


suggests  that  for  R  »  1,  where  the  inequality  is  strong  enough,  all  but  a 
negligible  fraction  of  the  total  dissipation  occurs  at  wavenuiribers  k  «  R  k  =  v  /v. 
We  shall  assume  that  this  is  so  and  verify  later  that  the  resulting  equations  do 
yield  a  dissipation  which  falls  off  rapidly  for  k  »  k,.  Then  for  k  in  the 
dissipation  as  well  as  the  inertial  range,  r(k,s)  ~  g(k, |s|)  is  given  by  (5«6), 
Let  us  anticipate,  on  the  basis  of  the  local  cascade  process  already  found  for 
the  inertial  range,  that  energy  transport  in  both  inertial  and  dissipation  ranges 
is  negligibly  dependent  on  direct  contributions  from  triad  interactions  with 
modes  outside  both  these  ranges.  Then  the  energy  balance  equation  (U.2)  takes 
the  asymptotic  form 


(7.1)   2vk^(k)  =  kjf  [k^a(k,p,q)E(p)-p^(k,p,q)E(k)]E(q)©j(k,p,q)-^3>, 

where 

[^  JlC^^o^s)  J-l(2v^Ps)  J^{2v^q3) 
vjcs     v^ps     v^qs 

^o      O  O  0 


(7.2)         e^{k,p,q)   = 


As  before,  we  have  taken  F(k)  to  be  negligible  for  k  »  k  .  It  is  clear, 
physically,  that  the  solution  of  (7.1)  cannot  be  unique,  since  we  have  not 
specified  the  power  input  or  dissipation  rate,  (In  effect,  we  have  pushed  the 
source  of  energy  supply  off  to  zeitj  wavenumber.)  To  specify  a  unique  solution 
we  may  requijre 
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(7.3)  I  2vk^(k)dk  =  e. 


If  we  write 


(7.Ua)  E(k)  -  (6V^)^/V3/2^(kA^), 


which  may  be  expressed  in  the  dimensionless  form 


(T.Ub)        E(k)/V^v  -  (kA^)-^/2^(kA^), 


we  readily  find  that  (?•!)  reduces  to  the  universal  equation 

ff 


(7.5)    2u'f(u) 


A  I- 


u2a(u,v,w)(2p/2f(.).A(..v.w)f(u)  yr^'^tWz[l't\^  ' 


U  U/  TV 


Where 

/°"Jt(2s)  Jt(2iiis)  J,  (2ns) 

(7.6)  K(m,n)  =  J^  J~-  -i~ ^^^  ds. 

Equation  (7.3)  takes  the  form 

(7.7)  f  u^/2f(u)du  =  \   . 

Instead  of  imposing  (7.7),  we  could  require  that  f(0)  have  the  value  determined 
in  Section  6  from  (6.2).  The  two  conditions  are  equivalent  because  the  nonlinear 
interaction  is  conservative;  the  energy  transported  through  the  inertial  range 
is  equal  to  the  total  viscous  dissipation. 

The  singularities  in  the  integrand  of  (7.5)  at  v  -  0  and  w  »  0  cancel, 
as  may  be  verified  from  (3.8)  and  (3.9).  This,  again,  expresses  the  localness 
of  the  cascade  process,  according  to  which  the  total  contribution  to  the  integrand 
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from  the  infinitesimal  neighborhood  of  these  points  is  infinitesimal.  We  may 
anticipate  that  f (u)  is  a  well-behaved  function  in  the  entire  range  0  <  ti  <  ■^ , 

It  was  pointed  out  in  Section  U  that  the  emission  term  in  the  energy 
balance  equation  is  proportional  to  E(k)  so  that,  if  the  spectrum  level  in  the 
immediate  neighborhood  of  k  were  suddenly  raised  above  its  equilibrium  value, 
the  result  would  be  a  decrease  of  the  net  nonlinear  input  to  the  neighborhood 
and  a  return  to  eqiiilibrium.  This  behavior  suggests  that  it  should  be  possible 
to  solve  (7»5)  by  a  corresponding  correction  or  iteration  procedure  in  which 

one  obtains  sin  improved  trial  function  as  an  appropriate  linear  combination  of 

2 

an  initial  trial  function  2u  f (u)  and  the  function  generated  as  the  right  side  of 

(7,5)  upon  substituting  the  initial  function.  It  may  be  expected,  however,  that 
considerable  practical  difficiilty  will  arise  because  of  the  singularities  in  the 
integrand. 

It  is  of  interest  at  this  point  to  discuss  further  the  relative  roles 
plajred  in  (5.U)  by  local  and  by  energy  containing  modes.  We  have  noted  that  for 
k  «  R  k  the  behavior  of  r(k,T)  and  g(k,T)  as  functions  of  t  ig  determined  by 
interactions  involving  the  energy  containing  region.  However,  although  contri- 
butions involving  the  energy  containing  region  dominate  the  right  side  of  {^»h) 
away  from  the  origin,  they  cancel  out  in  the  immediate  neighborhood  of  x  =  0,  as 
indicated  by  (5.5),  In  this  neighborhood,  r(k,T)  goes  to  zero  and  the  term 
vk  r(k,T),  which  here  dominates  the  left  side,  is  balanced  by  contributions  on 
the  right  nearly  entirely  from  local  triads.  At  t  «  0,  this  corresponds  to  the 
localness  of  energy  transport.  As  the  inequalities  k  «  k  «  R.k  become 
stronger,  the  range  of  t  about  the  origin  in  which  wholly  local  contributions 
dominate  becomes  a  smaller  fraction  of  (v  k)"  . 


-  31  - 


7.2.  Behavior  for  k  »  k 

Let  us  suppose  that  R  is  sufficiently  large  that  there  are  values  of 
k,  very  large  compared  to  k,,  for  which  (7.1)  is  still  valid.  There  appear  to 
be  two  general  possibilities  for  the  energy  supply  to  such  wavenumbers.  Either 
they  are  powered  by  an  essentially  local  cascade  process,  or  they  draw  power 
principally  from  wavenuinbers,  the  order  of  k,,  in  the  ref^ion  where  the  spectrum 
begins  to  fall  below  the  inertial  range  form.  We  might  perhaps  anticipate  that 
the  first  possibility  is  the  actual  fact.  In  the  inertial  range  we  have  seen 
that  the  energy  transport  is  essentially  local.  Let  us  now  imagine  a  part  of  the 
spectrum  to  be  depressed  by  local  viscous  drain.  It  would  appear  that  the  net 
input  from  much  lower  wavenumbers  should  be  affected  relatively  little,  since  to 
start  with  the  local  region  was  already  very  far  fixsm  equipartition  with  these 
wavenumbers.  On  the  other  hand,  the  local  decrease  in  the  level  of  excitation 
with  respect  to  that  of  closer  wavenumbers  should  act  in  the  direction  of 
increasing  relatively  more  the  net  input  from  these  latter  wavenumbers.  .These 
considerations  actually  must  be  considerably  elaborated  because  of  the  non- 
linearity  of  the  system,  but  it  is  not  difficult  to  verify  that  the  transport 
does  remain  essentially  local  in  the  dissipation  range  and  to  find  the  quali- 
tative shape  of  the  spectrum  for  k  »  k,. 

First,  it  is  easy  to  see  that  no  law  of  the  form  E(k) oc  k  is  possible 
for  k  »  k,.  Consider,  to  start,  the  assumption  that  the  transport  is  essentially 
local  so  that  the  right  side  of  (7.1)  is  dominated  by  contributions  for  which 
q  and  p  are  the  order  of  k.  Then  in  the  contributing  regions  9j(k,p,q)  is  of 
order  (v  k)~  ,  Since  a  and  b  are  dimensionless,  one  finds,  counting  up  the 
wavenumber  factors,  that  (7»l)  corresponds  to  the  order-of -magnitude  relation 
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(7.8)  2vkS;(k)  ~  ▼^"V[E(k)]^. 

As  k  increases  indefinitely  this  cannot  be  satisfied  by  any  power  law  except  the 
absurd  choice  n  =  0,  On  the  other  hand,  if  we  assttme  that  most  of  the  power  is 
supplied  by  interactions  involving  modes  of  wavenumber  ~k,,  then  in  place  of 
(7»8)  we  have 


(7.9)  2vk^E(k)  ~  v^"^^kE(k^)E(k), 

where  leading  terms  have  been  retained  on  the  right.  The  expression  k,k  appears 

2 

instead  of  k  because  one  of  the  pair  p,  q  is  of  order  k,  and  the  integration 

ranges  of  both  p  and  q  are  restricted  to  intervals  of  width  ~k,.  Equation 

(7.9)  cannot  be  satisfied  by  any  power  law  at  all.  Thus,  both  (7.8)  and  (7.9) 
are  inconsistent  with  the  initial  assumption  of  a  power  law,  under  which  they 
were  derived. 

The  impossibility  of  a  power  law  arises  from  the  fact  that  E  appears 

linearly  on  the  left  of  (7.1)  and  bi.li nearly  on  the  right,  while  the  equation 

is  homogeneous  in  the  wavenumbers.  This  suggests  that  E(k)  contains  a  factor 
-ck>^, 
e      (where  c  is  a  nua»rical  factor),  since  that  form  gives  the  product  law 

-cpA,  -cqA,    -c(p+q)A, 
e    *^e      -  e         .An  investigation  of  (7.1)  for  k  »  k^  shows  that 

it  is,  in  fact,  satisfied  asymptotically  by  a  function  of  the  form 

2  -ckAj 

(7.10)  E(k)  =•  const  (kA^)  e    *^        (k  »  k^). 

The  dominant  contributions  to  the  power  input  come  from  wavenumber  pairs  p,  q 
forming  very  flat  triangles  with  k,  as  would  be  anticipated  from  the  effect  of 
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the  exponential  factors  in  depressing  the  value  of  E(p)E(q)  for  psdrs  not  having 

this  property.  The  integration  thus  is  effectively  confined  to  an  area  of  the 

p,  q  plane  very  much  smaller  than  k  ,  and  in  this  area  E(p)E(q)  »  [E(k)]  .  The 

factor  (kA^)  »   which  might  seem  surprising  in  view  of  the  homogeneity  of  (7«l) 

in  wavenumbers ,  arises  because  of  the  restriction  of  integration  domain;  its 

form  also  is  affected  by  the  behavior  of  a(k,p,q),  which  vanishes  for  an  exactly 

flat  triangle.  As  we  partially  anticipated  on  general  grounds,  the  energy 

transport  is  localj  this  appears  as  a  consequence  of  the  fact  that  the  fraction 

of  P>  q  pairs  satisfying  p  +  q  ~  k,  and  such  that  p  or  q  ~  k,,  goes  to  zero 

with  k,A. 
d 

The  form  of  (7»10)  confirms  the  assumption  that  the  dissipation  is 
effectively  confined  to  wavenumbers  the  order  of  k,«  It  suggests  also  that 
each  decade  range,  say,  of  wavenumbers  k  »  k,  dissipates  most  of  the  energy 
supplied  to  it,  passing  on  an  ever-decreasing  fraction  to  the  next  decade. 
Corresponding  to  this  fact,  the  emission  terms  in  (7»l)  are  very  small c ompared 
to  the  absorption  term  for  k  »  k,. 

An  important  characteristic  of  (7.10)  is  that  it  suggests  the  existence 
of  the  nraments  f     k"E(k)dk  for  arbitrarily  high  n.  In  the  x  space  representation, 
this  implies  the  existence  of  mean  square  velocity  derivatives  of  «tll  orders  - 
a  property  which  seems  required  by  physical  intuition.  To  establish  this 
conclusion  more  fully  it  is  necessarj'  to  investigate  the  behavior  of  E(k)  for 
extremely  high  wavenumbers  k  »  R  k  ,  Although  we  have  not  determined  the  form 
of  r(k,T)  in  this  range,  the  energy  balance  equation  is  nevertheless  tractable 
because  the  function  6(k,p,q)  appearing  in  the  dominant  absorption  term  is 
determined  asymptotically  by  g(k,s)  -  exp(-vk  s).  An  analysis  of  the  type 
described  above  shows  that 
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(7.11)  E(k)    ^     AspK const )k]  (k  »  R  k  ), 


o  o 


1/3 

This  resTilt  appears  to  be  valid  even  for  R  ~    1»     The  proportionality 

constants,  however,   are  not  independent  of  R  , 

7«3»     Skewness  and  flatness  factors  of  the  distribution  of  3  u^/3x- 

1/3 

For  R  »  1,  the  skewness  factor 
o 

(7.12)  S^  =  ((du^/3x3^)^>/<(aS^/ax^)^>^/2 

is  given  in  terms  of  spectrum  moments  by  the  asymptotically  valid  relation 
(Batchelor,  1953,  p.  168) 

(7.13)  S^  -  -|(30)^/2v[j    k2E(k)dk  J    k^:(k)dk. 
By  (7.U)  and  (7.7),  we  find 

(7.1U)  S^   =  -  ^(l5)^/\-^/^  f"u^/2f(u)du. 

17         0     ^ 

We  Can  also  evaluate  S^  by  expressing  it  in  terms  of  k  space  triple  moments  and 
expressing  these  moments  in  terms  of  U(k,T)  and  g(k,r)  by  the  method  employed  in 
Paper  I  to  obtain  (3.7).  This  same  procedure  can  be  applied  in  general  to  the 

determination  of  the  skewness  factor  S  of  the  distribution  of  8  u,/3x,  .  It  is 

1/3 

not  difficult  to  show  that  for  R  '"^  »  1  the  result  depends  asymptotically  only 

on  the  dissipation  range  structure  and  that  upon  substituting  (5.3),  (5.6),  and 
(7,U)  in  the  resulting  expressions  we  obtain  the  asymptotically  valid  result 

(7.15)  S  =  A  R  '^Z^, 

^  n    n  o    ' 

where  the  A  are  universal  constants. 
n 
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The  R  dependence  indicated  by  (7»lh)  and  (7.l5)  is  quite  weak,  and 
it  is  impossible  to  estimate,  without  further  analysis,  how  strong  the  inequality 
R  '-^  »  1  has  to  be  for  the  asymptotic  behavior  not  to  be  masked  by  other 
effects. 

The  application  of  the  theory  to  fourth-order  moments,  which  is 
outlined  in  Paper  I,  yields  an  expression  in  terms  of  U(k,T)  and  g(k,T)  for  the 
flatness  factor  of  the  distribution  of  d  u,/3x,  .  It  can  be  seen  from  this 

expression  that  the  flatness  factor  approaches  for  each  n  some  universal  value 

1/3 

asymptotically  independent  of  R  ,  for  large  R 

7,ll,     A  visualization  of  the  dynamics  in  terms  of  linear  systems 

An  interesting  paradox  arises  when  we  attempt  to  interprete  the 
dynamics  of  huge  Reynolds  number  turbulence,  as  given  by  the  present  theorj',   in 
terms  of  familiar  linear  systems.     Equation  (6.?)  implies  that,   if  the  driving 
forces  were  removed,  the  characteristic  decay  time  for   the  energy  containing 
range  of  the  turbulence,  due  to  the  power  output  to  higher  wavenumbers,  would  be 
the  order  of  (v  k  )"  j  that  is,   the  energy  range  is  approximately  critically 
damped  by  a  resistive  dynamical  coupling  to  higher  wavenumbers. 

Now  let  us  go  up  the  wavenvimber  spectrum  considering  each  decade,  say, 
a  sub-system.     As  we  go  through  the  inertial  range,   the   characteristic  time  of 
each  sub-system  is  of  order  (v  k)"  ,   the  energy  in  each  sub-system  is  of  order 
(cv  )  '  k     '    ,  and  the  power  output  to  the  next  higher  sub-system  is   ~ e.     The 
"damping  factor"  of  the  sub-system,   the  fraction  of  energy  passed  on  in  a 
characteristic   "period",   is  of  order 


(V  k)-V(r.  )l/Vl/2    =     (k  A)l/^ 


6,- 
O 
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Thus,  the  damping  due  to  dynamical  coupling  goes  steadily  down  as  we  rise 
through  the  inertial  range. 

When  we  get  into  the  dissipation  range  k  >  k,,  a  similar  argument 
shows  that  the  damping  due  to  coupling  with  the  next  higher  sub-system  continues 
to  go  down,  and  with  great  rapidity  since  the  coupling  to  higher  modes  becomes 
ineffective  and  transfers  a  rapidly  decreasing  fraction  of  the  energy  in  a  sub- 
system. In  this  range,  as  in  the  inertial  range,  the  damping  factor  associated 
with  direct  viscous  dissipation  increases  linearly  with  k  (as  vk  /v  k)j  however, 

it  does  not  become  critical  until  k  ~  R  k  ,  which  is  far  above  the  range  of 

0  0*  ^ 

significant  dissipation  for  huge  Reynolds  numbers.  Thus  we  have  a  situation  in 
which  the  energy  containing  range,  where  viscosity  does  not  act,  is  approximately 
critically  damped,  while  the  ran^je  in  which  viscosity  disposes  of  the  energy  is 
very  lightly  damped!  As  we  go  up  above  R  k  ,  the  damping  due  to  dynamical 
coupling  continues  to  decrease  rapidly  (it  has  long  since  become  extremely  small), 
while  the  direct  viscous  damping  becomes,  finally,  very  large  and  dominates  the 
mechanics  of  the  sub-systems. 

The  behavior  described  seems  strange,  but  there  does  not  appear  to  be 
anything  either  physically  or  mathematically  inconsistent  about  it.  It  arises 
from  the  nonlinearitjr  and  complication  of  the  system,  and  the  consequent  arti- 
ficiality of  dividing  it  into  sub-systems.  The  characteristic  frequencies  of 
the  sub-systems  are  not  entirely  internal  parameters,  but  are  strongly  affected 
by  the  energy  range  mixing,  as  discussed  before.  We  must  also  keep  in  mind  that 
each  sub-system  serves  as  a  kind  of  frequency  multiplier  for  the  energy  supplied 
to  it» 

The  asymptotic  identity  of  time  correlation  and  response  functions  for 
the  inertial  and  dissipation  ranges  indicated  by  the  present  theory  is  of 
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considerable  fundamental  interest.  It  has  been  shown  for  many  types  of  conserva- 
tive systems  that  time  correlation  and  averaged  response  functions  are  identical 
in  an  equilibrium  (detailed  balance)  ensemble.  In  particular,  this  is  true  for 
a  dissipationless  system  closely  related  to  turbulence  (Paper  I,  appendix).  In 
the  present  case,  we  have  found  this  same  behavior  for  modes  which  are  very  far 
from  a  state  of  detailed  balance  and  which  dissipate  most  of  the  energy  lost  by 
the  system.  It  seems  likely  that  this  is  connected  with  the  fact  that  for  these 
modes  the  effective  damping  due  to  both  dynamical  coupling  and  direct  dissipation 
is  very  small. 

8.  COMPARISON  WITH  EXPERUffiNT 

8,1,  The  effect  of  deviations  from  asymptotic  conditions 

The  analysis  of  the  present  paper  has  been  based  on  exact  isotropy 
and  stationarity,  and  we  have  indicated  solutions  of  the  equations  only  for 
extremely  high  Reynolds  numbers.  It  is  not  possible  to  make  a  very  meaningful 
comparison  with  experiment  without  some  idea  of  the  theoretical  corrections 
implied  by  the  unavoidable  experimental  deviations  from  asymptotic  conditions. 
Unfortunately,  not  very  much  can  be  said  on  this  matter  at  the  present  stage. 
At  the  least,  a  semiquantitative  estimate  of  the  simplest  deviations  must  wait 
upon  evaluation  of  the  function  f(k/k^)  which  characterizes  the  asymptotic 
solution  itself.  For  the  present  we  can  attempt  only  a  very  qualitative  discussion. 

First,  we  may  consider  the  relation  between  our  fictional  stationary 
isotropic  turbulence  and  freely  decaying  isotropic  turbulence,  which  possibly 
can  be  closely  approximated  in  the  laboratory.  For  R^  sufficiently  high  that 
the  asymptotic  solution  developed  in  Section  7  has  a  range  of  accurate  validity. 
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it  seems  clear  that  the  behavior  of  this  range  should  be  negligibly  affected  by 
decay,  provided  enough  time  has  elapsed  for  the  high  wavenmnber  structure  to 
have  been  fully  established.  This  follows  from  the  fact  that  the  decay  of  v  is 
very  slow  compared  to  characteristic  periods  of  the  range.  The  spectrum  and 
time  correlation  for  the  stationary  case  therefore  should  apply  to  the  free 
decay  case,  provided  v  and  e  are  measured  within  a  time  short  compared  to  the 
decay  time.  Since  the  instantaneous  rate  of  viscous  dissipation  equals  the 
instantaneous  rate  of  decrease  of  kinetic  energy,  it  should  not  matter  whether 
6  is  measured  as  -3(4v  )/9t  or  as  /  2vkTJ(k)dk,  if  accurately  isotropic  and 

t  O  Jo 

homogeneous  conditions  have  been  established. 

At  lower  Reynolds  numbers  the  equivalence  of  the  stationary  and  decay 
cases  no  longer  follows,  and  we  must  consider  also  that  (7.1)  no  longer  is 
necessarily  an  accurate  expression  of  the  stationary  theory.  After  the  numerical 
constants  are  evaluated,  (6,$)  may  be  of  value  in  the  qualitative  understanding 
of  lower  Reynolds  numbers* 

The  role  of  anisotropy  in  the  present  theory  is  quite  interesting. 
Let  us  consider  huge  Reynolds  number  homogeneous  axisymmetric  turbulence.  For 
this  case  the  response  and  time  correlation  functions  for  high  wavevectors  can 
be  determined  in  close  analogy  to  the  isotropic  case,  yielding  similar  Bessel 
function  expressions.  It  is  easy  to  see,  however,  that  the  effective  mixing 
velocity,  replacing  v  in  the  arguments,  depends  on  the  angle  between  k  and  the 
axis  of  symmetry  and  is  maximum  when  k  is  in  the  direction  of  maximum  rms  velocity. 
The  energy  transport  among  a  triad  k,  p,  q  in  the  inertial  or  dissipation  ranges 
has  a  consequent  dependence  on  the  orientations  of  these  wavevectors  relative  to 
the  axis  of  symmetry.  Thus,  despite  the  localness  of  energy  cascade,  it  does 


-  39  - 


not  appear  likely  that  the  present  theory  yields  asymptotic  isotropy  of  energy 
partition  in  the  inertial  and  dissipation  ranges.     The  geometry  of  the  transport 
terms  is  sufficiently  complicated  that  further  investigation  is  needed  to  determine 
whether  the  anisotropy  in  the  9  factors  actually  results  in  appreciable  energy 
anisotropy  and,   if  so,  whether  it  is  of  the  same  or  opposite  sign  as  the  anisotropy 
of  the  rms  velocity  components. 

As  we  go  to  wavenumbers  k  >  R  k  ,   the  transport  terms  begin  to  be 

2 
dominated  by  isotropic  factors  of  the  form  exp(-vk  s)  in  the  response  functions. 

In  consequence  of  the  localness  of  the  cascade  process,  we  should  then  expect  an 

approach  to  energy  isotixjpy  at  these  high  wavenuitibers,   regardless  of  possible 

anisotropy  lower  in  the  spectrum.     This  should  be  true  regardless  of  the  size  of 


%' 


The  presence  of  inhomogeneity  as  well  as  anisotropy  brings  into 


question  the  basic  hypotheses  of  the  present  theory.     If,   for  example,  we  con- 
sider turbulence  confined  within  a  region  not  large  compared  to  k       ,   the  mode 
density  will  not  be  high,   and  we  cannot  expect  the  weak  dependence  property  to 
hold  accurately.     However,   the  weak  dependence  theory  may  be  expected  to  give  an 
exact  description  of  turbulence  hon»geneous  in  only  one  or  two  directions  -  for 
example,  flow  in  an  infinitely  long  pipe  or  channel.     To  obtain  closed,  deter- 
minate covariance  equations,  weak  dependence  need  be  assumed  only  for  wavevectors 
differing  in  their  axial  components;     the  phase  correlations  in  transverse 
directions  may  be  retained  and  determined,   along  with  the  mean  velocity  profile. 

Without  further  investigation,  extreme  ceution  must  be  used  in 
attempting  to  interpret  shear  flow  experiments  on  the  basis  of  homogeneous 
theory.     Laboratory  "spectrum"  measurements  involve  a  kind  of  "local  Fourier 
analysis"  and  give  results  which  depend  on  the  position  in  the  inhomogeneous 
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flow.     The  Fourier  components  used  in  the  present  theory  are  not  local  quantities j 
they  i^fer  to  the  entire  flow.     Without  further  justification,  we  cannot  divide 
into  local  regions  an  inhoraogeneous  flow,   in  which  there  are  necessarily  phase 
correlations  between  pairs  even  of  high-lying  Fourier  modes,   and  attempt  to 
treat  each  i^egion  more  or  less  independently.     It  is  conceivable  that  the  mixing 
action  of  a  strong  shear  velocity  in  one  part  of  the  flow  may  influence  the  energy 
balance  in  high  wavenuniber  Fourier  components  which  affect  a  laboratory  "spectrum" 
measurement  made  in  a  different  part  of  the  flow. 

Apart  from  the  interpretation  of  the  present  theory,   the  concept  of 
local  Fourier  analysis  is  not  a  simple  one.     Terms  like   "eddies  of  wavenumber  k", 
sometimes  used  in  discussing  turbulence,  imply  a  vagueness  of  distinction 
betvreen  x  space  and  k  space  concepts  which  can  lead  to  serious  inconsistencies. 

We  shall  now  proceed  to  disregard  the  cautions  just  presented  and 
compare  the  asymptotic  stationary,  isotropic  theory  with  real  experiments, 
including  shear  flow  experiments.     It  is  hoped  that  the  comparison  may  lead  to 
insights  regarding  the  questions  which  have  been  raised,     A  further  motivation 
is  that  the  only  other  course  open  at  present  would  appear  to  be   to  ignore 
experiments  entirely, 

8,2.     Turbulent  pipe  flow 

The  highest  laboratory  values  of  R     (>  10"^ )    appear  to  have  been 
obtained  in  shear  flows.     An  experimental  situation  for  which  very  careful 
measurements  have  been  taken  is  the  fully  developed  flow  in  a  long  circular  pipe 
(Laufer,  195U),     Measurements  were  made,  at  several  stations  along  the  pipe 
radius,  of  the  one-dimensional  spectra  ?f^(k^),  <^^ik^) ,  0-(k^)  of  the  energy  in 
the  axial,  radial,  and  circumferential  velocity  components,  respectively. 
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The  resoilts  of  the  measurements  on  the  axis  are  shown  in  Figure  ^<,     The  one- 
dimensional  spectrum  of  the  total  kinetic  energy  ii^(k_  )  »  0.  (k.  )  +  0-(k- )  +0_(k.  ), 
for  two  stations,  is  shown  in  Figure  6.     In  most  cases  the  points  shown  on  the 
latter  plot  are  not  experimental  points  but  are  computed  from  the  curves  drawn 
by  Laufer  through  the  original  experimental  points  for  the  three  component- 
spectra.     In  both  Figures  S  and  6,  wire  length  corrections  have  been  applied 
only  to  0^(k.).     The  values  of  R     '    ,  k  ,   and  k,  shown  on  Figure  6  are  computed 

2       1~2       ~2       ~  2 
from  V       »  •r'(u-     +  Sp     +  S,  )  and  estimates  of  local  values  of  e.     The  dissipa- 
tion rate  cannot  be  determined  precisely  from  the  measurements,  but  the  computed 

1/3 

parsmeters  are  not  very  sensitive  to  its  value.     R  and  k,  probably  are  not 

in  error  by  more  than  10%    and  k     probably  by  not  more  than  20  z^.     The  accuracy 
of  the  relative  values  of  the  component  spectra  may  be  estimated  as  about  l5  -  20 7o 
(Laufer,  private  communication). 

It  can  be  seen  from  Figure  6  that  the  total  energy  spectrum  follows  a 

-3/2 

k_      '     law  quite  closely  for  a  substantial  fraction  of  the  range  between  k     and 

k,.     If  estimated  wire-length  corrections   (Laufer,  private  communication)   are 
d 

-3/2 
applied  to  0p(k^)   and  0_(k^),   the  range  of  adherence  to  the  k,      '     law  is  some- 
what extended.     Within  the  range  where  <jiik~)  exhibits  -3/2  law  behavior,   it  can 
be  seen  from  Figure  5  that  there  are  substantial  differences  among  the  slopes  of 
the  component  spectra,   the  deficiency  in  excitation  of  u^  and  u,  at  low  wave- 
numbers  changing  to  an  excess,  above  isotropic  relation  to  u^, at  high  wavenumbers. 

Superficifilly  at  least,  the  behavior  would  seem  not  inconsistent  with 
what  might  be  expected  on  the  basis  of  the  asymptotic  isotropic  results.     The 
total  energy  cascade  appears  to  proceed  according  to   the  isotropic  law.     As  the 
cascade  goes  on,  there  appears  to  be  a  drain  of  energy  from  u^  to  Ug  and  u.  so 
that  the  component  spectra  have  respectively  greater  and  lesser  slopes  than  the 
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Figure  5«     Energy-component  spectra  0^(k^),  0Ak-),  0'iiK)  on  the  a^d-s  of  fully- 
turbulent  pipe  flow.     The  dashed  line  gives  0^^^^  ^3  calculated  from  0,(lc,) 
using  isotropic  relations.     The  pipe  radius  is  a  =  12,33  csm,  and  the  Reynolds 
number  based  on  a  and  the  mean  velocity  on  axis  is  5  x  10  ,     After  Laufer  (195U). 
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total  spectrum.  The  eventual  nonisotropic  partition  of  energy  at  high  wavenuribers 
does  not  appear  too  surprising  in  view  of  the  discussion  in  Section  8,1,  We  may 
note  that  R  k  ,  the  wavenuitiber  at  which  the  response  functions  must  begin  to 
become  isotropic,  is  far  above  the  range  plotted. 

In  the  original  interpretation  of  the  measurements  (Laufer,  195U),  it 
was  noted  that  0, (k^ )  follows  a  k,  '   law  very  closely  over  a  substantial  range. 
This  is  the  inertial  range  law  predicted  by  the  Kolmogorov  theory  for  isotropic 
turbulence,  and  the  result  was  taken  as  indicating  consistency  with  that  theory. 
The  comparison  does  not  seem  very  satisfying,  however,  as  Laufer  noted.  In  the 
center  of  the  inertial  range,  (2(^(k,)  represents  only  about  20  -  25'%  of  the 
total  energy  associated  with  k_  ,  and  the  other  components,  which  contain  most 
of  the  energy,  do  not  follow  the  -5/3  law  at  all.  In  making  a  comparison  with 
an  isotropic  theory  it  would  appear  more  reasonable  to  consider  the  total 
energy  spectrum,  as  we  have  done  above.  For  this  spectrum,  the  -3/2  law  appears 
to  give  an  optimum  straight-line  fit,  and  the  -5/3  law  appears  to  lie  outside 
the  experimental  error.  It  may  also  be  noted  that  the  very  fact  of  the  observed 
anisotropy  at  very  high  wavenumbers,  which  is  not  surprising  on  the  basis  of  the 
present  theory,  seems  a  complete  mystery  on  the  Kolmogorov  theory. 

It  should  be  emphasized  that  the  discussion  above  is  merely  suggestive. 
The  cautions  expressed  in  Section  8.1.  must  be  kept  firmly  in  mind. 

Measurements  in  the  boundary  layer  on  a  flat  plate  at  values  of 

1/3 

R  '  similar  to  those  in  the  pipe  flow  experiment  have  been  reported  by  Klebanoff 

(1955).  In  the  outer  part  of  the  boundary  layer  a  behavior  qualitatively  similar 
to  that  described  above  was  foundj  the  slope  of  9^-i(k, )  was  steeper,  and  that  of 
02(k-.)  flatter,  than  -3/2.  Unfortunately,  measurements  of  0o(^ )  were  not  made. 
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It  is  of  interest  that  the  anisotropy  at  very  high  wavenumbera  found  in  the 
pipe  flow  measurements  did  not  appear  in  this  experiment. 

Betchov  (195?)  has  also  reported  measurements  at  closely  comparable 
R    ,  made  in  a  shear  flow  of  novel  character.  Only  0-1  (k,)  was  measured,  and 
it  resembled  closely  the  corresponding  function  for  the  pipe  flow.  No  meaningful 
comparison  with  isotropic  theory  seems  possible  in  the  absence  of  data  on 
^^(k^)  and  0^(k^). 

8.3.  Grid  turbulence 

Grid  produced  turbulence  appears  to  represent  at  present  the  closest 
laboratory  approach  to  homogeneity  and  isotropy.  Although  Rejoiolds  numbers  high 
enough  to  produce  an  inertial  range  have  not  been  described  in  the  literature, 
it  is  of  some  interest  to  examine  to  what  extent  the  dissipation  range  spectrum 
at  moderate  Reynolds  numbers  obeys  the  scaling  indicated  by  (7.Ub).  Figure  7 
shows  the  results  of  one-dimensional  spectirun  measurements  by  Stewart  and 
Tovmsend  (19^1)  as  plotted  by  these  authors  on  the  basis  of  scaling  indicated 
by  the  Kolmogorov  theoiy.  Similar  measurements  have  been  reported  by  Liepmann, 
Laufer,  and  Liepmann  (I95l).  Figure  8  shows  the  curves  drawn  by  Stewart  and 
Townsend  through  their  data  points,  rescaled  according  to  (7.Ub).  The  plots 
have  been  prepared  so  that  the  curves  for  R^^,  =  52^0  are  congruent  in  each  case. 
The  values  of  R^.  correspond  to  the  values  shown  in  the  table  below: 
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Figure  8.     The  curves  in  Figure  7  rescaled  according  to  the  present  theory 
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It  will  be  noticed  that  the  curves  of  both  k.  0n (k^)  and  k^^0^(k^) 

fall  substantially  closer  together  with  the  new  scaling  than  with  the  old.  In 
particular,  the  systematic  shift  with  TL.  of  the  horizontal  position  of  the  maxLiBa 

seems  largely  absent  with  the  new  scaling*  It  should  be  noted  that,  while  R  ^ 

2 

has  a  definite  theoretical  significance  as  the  wavenuniber  at  which  vk  =  v  k, 

k  is  only  nominally  the  wavenuiriber  characterizing  the  energy  containing  range. 

In  this  experiment  most  of  the  energy  lies  below  k^  in  each  case.  From  the 

values  given  in  the  table  it  would  seem  not  implausible  that,  particularly  in 

the  case  of  k-ViC^)*  "tbe  functions  should  bear  some  relation  to  the  asymptotic 

forms  for  a  substantial  part  of  the  k  range  plotted.  However,  the  values  of 

R  '^  are  too  low,  and  the  experimental  scatter  too  high,  to  permit  attaching 
o 

very  much  significance  to  the  comparison.  It  must  also  be  noted  that  detailed 
measurements  on  the  isotropy  of  the  flow  were  not  presented* 


9,  COMPARISON  WITH  OTHER  THBDRIES 

9,1,  The  Kolmogorov  theory 

The  Kolmogorov  theorj^  (Kolmogorov,  19lllj  Onsager,  19hSi   v,  Weizsacker, 
I9U85  reviewed  by  Batchelor,  19^3,  chap,  6)  has  occupied  a  central  place  in 
thinking  on  turbulence  in  the  past  decade  because  of  the  intuitive  appeal  of  its 
assumptions,  the  economy  of  its  methods,  and  the  approximate  empirical  support  for 
its  predictions.  Certain  predictions  of  the  present  weak  dependence  theory  appear 
to  differ  only  slightly  from  those  of  the  Kolmogorov  theory:  the  spectrum  law 
in  the  inertial  range  by  (kA^)    and  the  characteristic  scaling  wavenuniber  in 
the  dissipation  range  by  R  *"  '   ,  A  detailed  comparison  of  the  two  theories  is 
not  possible  because  one  is  wholly  quantitative  and  the  other  consists  of  dimen- 
sional arguments.  However,  we  shall  see  that  their  postulatorj""  bases  are  funda- 
mentally in  conflict. 
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Kolmoeprov's  original  assertion  was  essentially  that  the  internal 
dynamics  of  the  sufficiently  fine-scale  structure  in  x  space,  at  high  Reynolds 
number,  should  be  independent  of  the  large  scale  motionj  the  latter  should,  in 
effect,  merely  convect  local  regions  bodily.  The  restatement  of  this  hypothesis 
by  V.  Weissacker  holds  that  the  dynamical  interaction  of  the  Fourier  modes  u(k) 
is  effectively  local  in  wavenomber  space  and  that  low  k  modes  exert  on  sufficiently 
high  k  modes  only  an  action  effectively  equivalent  to  convection  by  a  suitably 
random  unifom  velocity  field.  This  implies  that  the  transport  of  energy  from 
low  to  high  wavenunibers  proceeds  by  a  cascade  process  and  that  the  internal 
mechanism  of  the  cascade  process  at  high  k  is  quite  independent  of  the  original 
source  of  the  energy.  Then  it  can  be  seen  that  the  spectrum  in  the  inertial 
range  depends  only  on  e,  whence,  by  dimensional  analysis,  E(k)  «=  e  '  k    .  In 
contrast,  our  inertial  range  spectrum  (6.3)  involves  the  energy  range  parameter 

^o- 

It  is  not  hard  to  see  that  the  weak  dependence  hypothesis  introduced 

in  Section  2  is  inconsistent  with  the  idea  that  the  low  k  motion  has  the  con- 
vective  effect  of  a  random  uniform  velocity  field  on  the  high  k  dynamics.     As  an 

illustration,   consider  an  ensemble  of  one-dimensional  flows  consisting  of  a 

2 

Gaussianly  distributed  uniform  velocity  v,  with  mean  square  v     ,   and  a  superposed 

zero-mean  fluctuation  field  u(x)  =  Z  u(k)e       ,  u(-k)   =  u^(k).     For  simplicity, 

k 

let  us  suppose  that  the  u(k)  are  sufficiently  weak  that  nonlinear  interaction 
among  them  is  negligible.  Let  us  assume  that  at  t  =  0  the  various  u(k,0)  are 
statistically  independent  of  each  other  and  of  the  uniform  field  v.  Then  for 
t  >  0  we  have 

(9.1)  u(k,t)  =  u(k,0)e~^''^, 
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which  expresses  the  conrective  action  of  v.  If  we  now  imagine  that  at  time  t' 
some  external  agency  produced  a  small  impulsive  disturbance  in  u(k),  this 
disturbance  would  simply  remain  and  be  convected  also.  Therefore  the  impulse 
response  function  for  u(k)  is  simply 

(9.2)  ^(k5t,t')  -  e"^""^^^"*'^     (t  >  f). 

The  normalized  time  correlation  function  and  averaged  impulse  response  function 
for  u(k)  are  then 


(9^3)         <u(k,t)u^kt.r)>^  ^(^  ^j  „  <^ivkr)  ^  ^^_  1^2k^2), 


(9.1)        <i^(kjt+T,t)>=  g(k,T)  -  exp(-|v^VT^)     (^>0), 

where  we  have  used  the  assumption  that  the  v  distribution  is  Gaussian, 
Now  let  us  consider  the  moment 


M  ^<^(k;t+ s,t)u(p,t)u(q,t)u  (p,t-t-3)u  (q,t-<-s)> 
(Mp,0)l^><h(q,0)|^> 

By  direct  evaluation  from  (9.1)  and  (9.2)  we  have 

(9.5)  M  =  exp  [-  |v  ^(k-q-p)^s^j  (s  >  O), 

where  the  independence  of  the  u(k,0)  has  been  used.     If,  instead,  the  weak 
dejDendence  hjrpothesis  were  asserted,  we  would  have  from  (9.3)   and  (9.h) 

Note  that  the  ensemble  averages  here  are  not  equivalent  to  neighborhood 
averages.     This  represents  a  fundamental  difficulty  in  the  explicit  formulation 
of  the  Kolmogorov-v,  Weizsacker  concepts. 
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(9.6)  M  ^   g(k,3)r(p,3)r(q,3)  =  exp[- ^v^^(k^  +  p^ +q2)3^J    (3  >  O). 

For  k-p-q  =»  0,  which  corresponds  to  the  triangle  relation  of  wavevectors  in  our 
three  dimensional  theory,  the  weak  dependence  hjrpothesis  would  give  a  very  wrong 
answer  indeed.  It  is  clear  that,  in  the  case  taken,  statistical  independence  of 
simultaneous  quantities  does  not  imply  weak  dependence  of  ron-siimiltaneous 
quantities. 

The  analysis  just  given  indicates  why,  in  actual  three-dimensional 
turbulence,  the  weak  dependence  theory  gives  a  dependence  of  high  wavemunber 
dynamics  on  v  ,  in  conflict  with  the  Kolmogorov  hypothesis.  The  9  factors  in 
the  energy  balance  equation  (7.1)  are  derived  from  moments  of  the  type  just 
considered.  For  p  +  q  =  k  they  could  not  have  a  dependence  on  v  ,  as  indicated 
by  (7.2),  if  the  energy  containing  modes  effectively  exerted  on  high  wavenumbers 
a  simple  convective  action  as  implied  by  the  Kolmogorov  theory.  It  is  of  some 
interest,  in  this  connection,  arbitrarily  to  replace  the  d  factors  in  (7.1)  by 
expressions  of  the  form 

.00 

(9.7)  9j'(k,p,s)  =  J  r(£aj^s)r(co  s)r(m  s)ds, 

where 


oi^  =  k[kE(k)]^/2 


is  a   "local"  characteristic  mode  frequency  as  suggested  by  the  Kolmogorov  theory. 

2/3  -5/3 
With  this  change,   (7.1)  leads  to  an  inertial  range  law  of  form  E(k)  oc    e  '   k     '   , 

and  we  find  a^cfz  t-'^\^'^ » 
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The  action  of  low  k  modes  on  high  k  modes  implied  by  the  weak  dependence 
theory  can  be  described  as  follows.  Let  us  suppose  that  at  t  »  0  the  various 
u(k)  are  weakly  interdependent,  and  let  us  then  switch  on  the  dynamical  inter- 
action. We  notice  fi^m  (3.1)  that  there  is  an  essential  difference  between  the 
actions  of  modes  of  low,  no n- zero  wavenumber  and  modes  of  zero  wavenuniber  (uniform 
field).  A  mode  k'  »  0  couples  mode  k  to  itself,  while  a  mode  k^'  ^  0  (but  however 
small)  couples  mode  k  to  another  mode  with  which  it  has  only  a  weak  statistical 
relation.  Now  let  us  take  three  high  wavenuniber  modes  k,  p,  q  and  consider  the 
effect  upon  them  of  the  presence  of  a  large  number  of  excited  modes  of  non-zero 
wavenumbers  «  k,  p,  q.  The  three  modes  are  not  coupled  each  to  itself,  yielding 
coherent  time  dependencies  of  the  tj-pe  (9.1).  Instead,  each  mode  has  mixed  into 
it  energy  from  many  closely  adjacent  modes  with  which  it  has  only  weak  statistical 
relations  and  which,  in  turn,  are  only  weakly  related  to  the  modes  being  mixed 
into  the  other  two  modes.  In  this  situation,  the  cancellation  of  time  factors 
expressed  by  (9.5)  does  not  occur.  It  can  be  shown  that  many  low  k  modes  must 
contribute  to  the  mixing  action  if  the  picture  is  to  be  self-consistent.  This, 
of  course,  is  ensured  by  our  limit  L-^oo, 

It  should  be  noted  that  the  presence  of  a  true  uniform  (zero -wave number) 

velocity  field  does  not  invalidate  the  weak  dependence  hypothesis.  A  uniform 

iv-k 
field  V  simply  adds  to  all  the  correlation  and  response  functions  a  factor  e  

or  e ,  as  is  clear  from  (3»l).  Averaging  over  an  ensemble  of  true  uniform 

fields  3:3  inadmissible. 

The  difference  of  action  between  uniform  fields  and  low  k  fields  may 

be  described  in  x  space  as  follows.  Suppose  our  large  box  of  dimension  L  is 

filled  with  initially  weakly  dependent  high  k  excitation.  If  a  low  k  motion 

exists,  it  will  shear  and  distort  each  initial  harmonic  component|  this  is  just 
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the  mixing  action  described  above.     The  shearing  distortion  appears  locally  to 
be  very  small,  but  it  is  clear  from  (3«l),  snd  the  discussion  above,   that  the 
mix3-n2  of  adjacent  high  k  modes  proceeds  at  a  rate  which  depends  on  the  magnitude 
of  the  velocity  and  not  of  the  shear  associated  with  the  low  k  motion.     The 
locally  tiny  shearing  actions  add  up  to  an  important  cumulative  effect.     This 
points  up  the  fact,   noted  earlier,   that  it  is  necessary  to  use  great  care  in 
combining  concepts   from  the  x  space  and  k  space  i*epresentations.     Quite  apart 
from  whether  it  is  profitable  to   think  of  turbulence  as  consisting  of  local 
regions  in  x  space,   the  use  of  Fourier  components  gives  a  representation  of  the 
problem  which  is  intiinsically  nonlocal  in  x  space.     Each  u(k)  refers   to  the 
entire  box. 

Turning  from  questions  of  representation  to  physics,  we  should  note 
at  this  point  that  the  local  properties  of  fine-scale  dynamics  appealed  to  in 
the  original  Kolmogorov  hypothesis  hardly  seem  to  be  implied  by  the  qualitative 
empirical  behavior  of  turbulent  fluids.     The  fine  structure  at  high  Reynolds 
numbers  does  not  consist  of  a  uniform  distribution  of  little  random  blobs  but 
rather  of  a  substantially  intermittent  distribution  of  sharply  defined,  extended, 
and  tangled  vortex  filaments  and  sheets   (Batchelor,  19^3*  Section  8.U).     The 
transfer  of  energy  from  low  to  high  wavenurabers  appears  to  involve,  in  x  space, 
the  stretching  of  these  filaments  and  sheets.     This  behavior  is  readily  observed 
by  stirring  ink  into  a  bathtub  of  watef  (R     >  10^  is  easily  achieved),     A  typical 
filament  appears  to  extend  throughout  a  substantial  part  of  the   total  turbulent 
domain. 

The  spatial  structure  described  does  not  seem  consistent  with  the 
idea  of  fine  scale  isotropy  within  local  regions  asserted  by  Kolmogorov  (19U1). 
Isotropy,  when  it  exists,   appears   to  be  a  statistical  property  shown  by  averages 
over  large  space  domains  containing  many  local  regions.     The  fact  that  the  space 
correlation  may  be  isotropic  for  small  separations  does  not  at  all  imply  isotropy 
within  local  resaons.     The  confusion  of  these  two  properties  leads  to  unwarranted 
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dynamical  assumptions  and  has  constituted  a  fundamental  illogic  in  the  use  of 
the  Kolmogorov  ideas.     In  general,  isotropy  for  small  correlation  distances 
implies  isotropy  of  local  averages  in  k  space  (see  Section  2),  but  only  because 
the  u(k)  give  a  highly  nonlocal  representation  of  the  x  space  structure. 

It  is  of  interest  to  note,   finally,   that  the  preservation  of  incon^jres- 
sibility  by  the  pressure  forces  during  the  stretching  of  filaments  implies 
analytically  certain  long-range  ordering  effects  in  the  fluid  (Batchelor  and 
Proudman,  1956). 

9,2.     Heiseriberg's  heuristic  theory 

Several  authors  have  proposed  models  for  the  transport  power   TKk) 
based  on  simple  physical  analogies   (Batchelor,  1953 >   chap,  6),     Of  these,   the 
heuristic  theory  of  Heisenberg  (19U8)  has  proved  the  most  popular.     Heisenberg 
assumed  that    IT(k)   could  be  represented  as  a  power  dissipation  by  all  wavenumbers 
<  k  due  to  the  action  of  an  effective   "eddy  viscosity*  which  represented  their 
coupling  to  all  wavenumbers  >  k.     In  analogy  to   the  expression  for  dissipation 
by  real  viscosity,  he  took 

(9.8)  n(k)-|    2|kJ    [q"^(q)]^/%lp^(p)dp, 

where  le  is  a  disposable  numerical  constant.     The  eddy  viscosity  |(     was  obtained 
by  dimensional  arguments  on  the  basis  of  the  Kolmogorov-v,  Weizaacker  hypothesis 

of  the  localness  of  k  space  dynamics.     Because  of  its  dimensional  structure,   (9.8) 
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leads  to  a  k     '     inertial  range  spectrum  law.     At  very  high  k  it  leads   to  a  k 

law,  implying  the  non-existence  of  mean-square  third-  and  higher-order  velocity 

derivatives.     The  possible  relation  of  (9,8)  to   (2,2)  has  been  investigated  by 

Tchen  (195U). 
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The  comparison  of  the  eddy  viscosity  and  weak  dependence  theories  is 
made  more  illuminating  by  first  modifying  (9.8)   to  yield  a  k     '      inertial  range 
lawj  Heisenberg's  basic  physical  assximption  seems  independent  of  the  conflict 
between  the  weak  dependence  and  Kolmogorov  theories.     Accordingly  we  shall 
consider  the  transport  function 

(9.9)  n(k)  -  J^    ^y\'^j    q"^(q)dqlp^(p)dp, 

where  k'  is  a  new  disposable  numerical  constant.  The  eddy  viscosity  j  f  is  here 
obtained  by  dimensional  considerations  under  the  requirement  that  it  depend 
inversely  on  v  ,  as  implied  by  our  previous  considerations.  The  rationality  of 

(9.9)  makes  it  perhaps  more  satisfying  intuitively  than  (9.8). 

From  (9.9)  we  immediately  obtain 

(9.10)  -  ^^   =  2Kk(v^k)-^[['  i.  E(k)E(p)  f   -  I  VE(k)E(p)  ^  j  . 

This  equation  has  been  written  in  a  form  to  facilitate  comparison  with  -dri(k)/3k 
as  given  by  the  right  side  of  (7.1).  The  time  factor  (v  k)~  corresponds  to  the 
6,  factor  in  (7.1).  The  most  immediate  qualitative  difference  between  the  two 
expressions  is  that  (9.10)  expresses  the  total  interaction  as  an  integral  over 
pair  interactions  (k,p)  rather  than  the  triad  interactions  (k,p,q)  which  enter 
(7.1) i  and  vfhich  are  fundamental  to  the  structure  of  the  original  equation  of 
motion  (3.1).  A  further  difference  concerns  the  sharp  distinction  in  (9.10) 
between  the  interactions  with  modes  greater  than  or  less  than  k.  The  absorption 
term  involves  only  p  <  k  and  the  emission  term  only  p  >  k.  In  (7.l)>  both  sorts 
of  modes  enter  both  terms,  and  the  net  input  appears  as  the  balance  of  both 
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ingoing  and  outgoing  contributions  from  every  triad  interaction}  this  more 
closely  resembles  other  dissipative  equilibrium  situations  in  statistical  physics. 

E(k)  appears  as  a  factor  in  both  the  absorption  and  emission  terms  of 
(9.3fi)j  suggesting  that  this  transport  expression  does  not  give  the  tendency  to 
return  to  equilibrium,  upon  a  perturbation  of  the  spectrum  in  the  neighborhood 
of  k,  that  was  noted  in  the  discussion  of  (7.1).  This  peculiarity  shows  up  more 
clearly  in  (9«9)j  and  also  in  the  unmodified  form  (9.8),  If  the  equilibrium  is 
disturbed  by  an  initial  increase  of  spectrum  level  in  the  range  >  k,  then  by 
(9.8)  or  {9*9)   the  effective  eddy  viscosity  acting  on  wavenumbers  <  k  is  in- 
creased. This  will  lead  to  an  increased  flow  of  energy  from  wavenumbers  <  k  to 
those  >  k  which  will  tend  to  increase  further  the  eddy  viscosity  acting  on  the 
lower  wavenunibers .  A  possible  instability  is  indicated.  On  general  statistical 
mechanical  grounds,  we  should  expect  precisely  contrary  behaviorj  an  increase  of 
excitation  in  modes  >  k  should  lead  to  a  decreased  flow  of  energy  to  them  from 
modes  <  k,  as  on  the  weak  dependence  theoiy. 

It  can  be  seen  from  (9.10)  that  for  k  very  large  compared  to  the 
wavenumbers  which  contribute  appreciably  to  the  dissipation  integral 
e  =  2v  I  p  £l(p)dp,  the  absorption  teirc  is  dominated  by  contributions  p  «  k  in 
the  principal  (real)  viscous  dissipation  region.  Thus,  the  transport  does  not 
proceed  by  local  cascade  for  such  k,  as  it  does  in  the  weak  dependence  theory. 

Most  of  the  comparisons  made  above  can  also  be  made  on  the  basis  of 
(9.8),  but  with  less  clarity.  For  this  case,  the  irrationality  of  the  eddy 
viscosity  expression  results  in  greater  asymmetry  between  absorption  and  emission 
terms. 
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9.3  •     The  analsrbical  theories  of  Helseiiberg  and  Chandrasekhar 

The  hypothesis  that  fourth-order  niDments  of  the  distribution  of  the 

u(k)   are  related  to  second-order  moments  as  in  an  independent  distribution 

implies  normal  relation  of  the  corresponding  x  space  moments.     It  is  a  much 

stronger  statement  than,   and  contains,   the  hypothesis  of  weak  dependence  for 

these  moments.     Fix)m  this  it  is  clear  that  the  independent  relation  hypothesis 

is  also  inconsistent  with  the  Kolmogorov  theory.     As  discussed  in  Section  2,   the 

independent  relation  hypothesis  is  inadmissible  on  both  empirical  and  theoretical 

grounds.     It  nevertheless  is  of  interest  to  compare     the  analytical  theories  of 

Heisenberg  (19U8,  Section  5)   and  Chandrasekhar  (1955),  which  are  founded  upon  it, 

with  the  present  theorj"". 

Heisenberg 's  analytical  theorj''  properly  is  appropriate  to  freely 

decaying  turbulence,  but  its  essential  structure  may  be  discussed  in  terns  of 

our  present  stationary  specialization  of  formalism.      (For  discussion  of  the 

inertial  and  dissipation  ranges  at  high  Reynolds  number,   the  stationary  and 

decay  cases  are  equivalent.)     Tlie  theory  involves  first  expressing  the  triple 

moment  S(k,'c),  defined  by  (3.5),   in  terms  of  fourth-order  moments  by  substituting 

for  u  (q.t)  the  bilinear  expression  obtained  by  integrating  (2,2)  with  respect 
m  — 

to  time.     Next,   the  fourth-order  moments  are  expressed  in  terms  of  U(k,Tr)  by 
means  of  the  independent  relation  hypothesis.     This  results  in  a  closed  equation 
for  U(k,T).     The  analysis  was  carried  out  for  v  =  0,  but  can  be  straightforwardly 
generalized  to  v  ^0. 

The  expression  found  for  S(k,'r)   can  be  obtained  from  our  result  (3.7) 

2  2 

by  1)  replacing  the  response  functions  g(p,s),  g(q,s)  with  exp(-vp  s),  exp(-vq  s)j 

2)  discarding  the  terms  involving  a(k,p,q)j  3)  retaining  from  the  terms  involving 
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b(k,p,q)  only  a  particular  linear  combination  of  parts  which  are  antisymmetric 
in  k  and  p  for  t  =  0.  The  alteration  of  the  response  functions  is  equivalent  to 
ignoring  the  effect  on  them  of  the  nonlinear  interaction.  The  discarding  of 
terms  amounts,  on  the  basis  of  the  weak  dependence  theory,  to  throwing  out  most 
of  the  interaction  altogether.  This  shows  up  most  strikingly  if  one  evaluates 
the  contribution  to  r(k,T),  for  high  k,  due  to  energy  range  mixing,  after  the 
fashion  in  which  (5.5)  and  (5.6)  were  obtained.  One  finds  that  this  contribution 
vanishes  as  k/k  -^  °^  .  This  result  apparentlj'-  led  Heisenberg  to  a  curious  chain 
of  reasoning.  Independence  of  r(k,-r)  on  v  clearly  is  absurd  in  a  stationary 
coordinate  system,  but  it  is  just  the  behavior  one  might  expect,  on  the  basis  of 
the  Kolmogorov  -  v,  Weizsacker  considerations,  in  a  sojrt  of  quasi-Lagrangian 
cooixiinate  sj-^tem  moving  with  the  low  k  motion.  Heisenberg  concluded,  therefore, 
that  possibly  the  independent  relation  hypothesis,  and  consequently  his  statistical 
equation  of  motion  (Heisenberg,  19U8,  equation  (83)),  might  really  be  appropi-iate 
to  such  a  coordinate  system.  He  proceeded,  on  this  bssis,  to  seek  solutions 
r(k,T)  of  a  form  similar  to  that  used  in  (9.7).  There  seems  to  be  no  analytical 
support  whatever  for  this. 

The  theory  of  Chandrasekhar  (1955)  wss  developed  by  analysis  in  x 
space.  Expressed  in  k  space,  it  differs  from  that  of  Heisenberg,  as  generalized 
above,  only  in  that  S(k,T)  is  expressed  in  terms  of  fourth-order  moments  by 
substituting  a  bilinear  integral  expression  for  u.  (k,t-T;)  rsther  than  for 
u  (q,t).  This  gives,  however,  quite  a  different  result  from  Heisenberg's, 
thereby  indicating  the  intrinsic  inconsistency  of  the  independent  relation 
hypothesis,  Chandrasekhar's  result  may  be  obtained  from  (3.7)  by  replacing 
g(k,s)  with  exp(-vk  s)  and  discarding  all  the  contributions  involving  b(k,p,q). 
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The  terras  retained  are  all  rigorously  positive-definite  for  t  =  0,  and  the 
theory  consequently  does  not  yield  conservation  of  energy  by  the  interaction 
(Kraichnan,  195?)  as  does  Heisenberg's  prescription.     However,   the  theory  has 
the  mathematical  virtue  of  being  expressible  by  a  partial  differential  rather 
than  an  integro-differential  equation  in  x  space. 

The  inconsistency  of  the  independent  relation  hypothesis  with  the 

Kolmogorov  theory  is  clearly  illustrated  by  the  attempt  of  Chandrasekher  (1956) 

-5/3 

to  reconcile  his  theory  with  the  k     '     inertial  range  law.     In  order  to  obtain 

the  formal  possibility  of  a  space  correlation  function  consistent  with  this  law, 
it  was  found  necessary  to  assume  that  the  full  space-time  correlation  function 
for  small  space-time  separations  was  independent  of  v  ,     This  conflicts  seriously 
with  the  Kolmogorov  theory,  which,  as  we  noted,  reouires  such  behavior  for  the 
correlation  function  in  a  quasi-Lagrangian,   not  a Eulerian  coordinate  system.     It 
also  conflicts  with  the  basic  physics  of  the  problem. 

Our  discussion  has  indicated  that  in  both  Heisenberg's  and 
Chandrasekhar's  theories  the  use  of  the  independent  relation  hypothesis  throws 
away  substantial  parts  of  the  statistical  interdependence  and  dynamical  coupling 
of  the  Fourier  modes.     It  should  be  noted  that  essentially  equally  significant 
omissions  could  be  expected  of  possible  theories  made  determinate  by  assuming 
that  certain  sixth,  eighth,  or  higher  order  cxunulants  vanish  instead  of  fourth 
order.     As  the  order  increases,   the  number  of  the  neglected  individual  cross- 
moment  cumulants  entering  physically  significant  expressions  rises  rapidly,  and 
their  total  contribution  cannot  be  ejqjected  to  decrease  with  increase  of  order. 
This  situation  is  connected  with  the  nonlinearity  of  (3,1) >   through  which  phase 
relations  of  large  groups  of  modes  are  established  by  networks  of  triad 
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interactions.  The  weak  dependence  theory  retains  cumulants  of  all  orders.  Tbe 
presence  of  g  in  the  expressions  for  third,  fourth,  and  higher  order  moments 
(Paper  I)  indicates  the  effect,  on  individual  cross  moments,  of  the  inter-relation 
of  all  the  modes. 

9,14..  The  theories  of  Proudman  and  Reid,  and  Tatsumi 

Proudman  and  Reid  (I951i)  and  Tatstimi  (195?)  assert  the  independent 
relation  hypothesis  only  for  the  distribution  of  simultaneous  values  of  the 
u(k,t).  They  are  able,  thereby,  to  obtain  closed  equations  for  the  time  change 
of  E(k)  in  freely  decaying  isotropic  turbulence.  No  predictions  are  obtained 
concerning  time  correlations.  The  analytical  structure  of  these  theories  is 
complicated,  and  we  shall  only  indicate  very  briefly  the  relation  with  the  present 
theorj'  for  the  inertial  range  at  very  high  Rejmolds  numbers.  It  should  be  noted 
that  tliese  theories  were  primarily  intended  to  describe  the  energy  containing 
range. 

We  shall  assume  that  there  exists  an  inertial  range,  far  above  the 
energy  containing  range,  in  which  energy  transport  takes  place  by  local  cascade 
and  in  which  the  direct  effects  of  viscosity  are  negligible.  In  this  range  we 
shall  assume  that  the  spectrum  obeys  a  k"  law.  Then  it  is  not  difficult  to 
verify  that  the  hypothesis  of  Proudman  and  Reid,  and  Tatsumi  leads  to  an  asymptotic 
expression  for  -3n(k)/8k  at  time  t  which  differs  from  the  right  side  of  (7.1) 
only  in  that  e_(k,p,q)  is  replaced  by 
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't 


E(k,3) 
E(k,t} 


ds. 


where  E(k,t)  is  the  value  of  E(k)   at  time  t, 
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From  empirical  knowledge  we  know  that  t  must  be  of  order  (v  k  )"  , 

O  0    0' 

the  overall  decay  time  of  the  turbulence,  whereas  we  have  seen  that  9^  is  of 
order  (v  k)"  for  significant  triad  interactions.  If  arguments  are  followed 
similar  to  those  which  led  to  and  verified  (6,3),  we  now  find,   instead, 


(9.11)  E(k)  ~    {cvk  )^/V2 

o   o 

for  the  asymptotic  inertial  range  spectinm.  This  result  involves  a  contradiction; 

2 

the  integral  of  2vk  E(k)  over  the  inertial  range  now  diverges  toward  low  wave- 
numbers,  a  fact  which  is  inconsistent  with  the  initial  assumption  that  viscosity 
effects  were  negligible.     This  suggests  that  the  theories  of  Proudman  and  Reid, 
and  Tatsumi  do  not  yield  an  inertial  range,     A  more  detailed  investigation 
confirms  this  conclusion. 

The  result  found  is  not  too  surprising.     The  independent  relation 
hypothesis  is  known  empirically  to  be  invalid  for  high  wavenumbers  (Batchelor, 
19^3,  Section  8.2).     The  appearance  of  t     instead  of  9j  in  the  transport 
expression  would  seem  to  indicate  that  the  effect  of  ignoring  the  existence  of 
fourth-order  cumulants  is  to  discard  most  of  the  effect  of  the  dynamical  inter- 
action on  the  time  correlation  and  response  functions  which  enter  9. 
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